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FOREWORD 
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34686 entitled "Bearing Tester Data Compilation, Analysis, and Reporting and Bearing 
Math Modeling" for the George C. Marshall Space Flight Center of the National Aero- 
nautics and Space Administration. The work was administered under the technical 
direction of the Materials and Processes Laboratory, Engineering Physics Division of 
the George C. Marshall Space Flight Center with Mr. Fred J. Dolan acting as project 
manager. 

This report describes the work performed by Spectra Research Systems Southeast- 
ern Operations during the July 1982 - April 1984 period. Mr. Rodney Bradford was the 
Spectra Research Systems Project Manager. The project technical staff was directed 
by Mr. Joseph C. Cody who was the Project Engineer. A listing of the key project 
staff members is shown below: 
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Ms. Linda S. New 
Mr. M. Anthony Stone 
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1.0 INTRODUCTION 


This report describes the work accomplished during the contract period July 1982 
- April 1984 in support of the MSFC Bearing and Materials Tester (BMT) Program. The 
objective of this activity is to support the development and operation of the BMT 
including data reduction and evaluation. Since the Shuttle Main Engine (SJ4E) 
turbopump bearings operate in an environment considerably more severe than conven- 
tional bearing systems, traditional analysis methods and bearing life models are no 
directly applicable. 

The MSFC BMT Program is therefore directed toward a better understanding of the 
various parameters that effect/or determine the SSME turbopump bearing operational 
characteristics and service life and to develop and verify design tool s applicable to 
these systems. Support of this program involves a broad spectrum of engineering 
analysis activities including static analysis of the BMT Shaft Bearing System, high- 
speed bearing system analysis, contact stress evaluation, bearing failure mode eval- 
uation and thermal modeling of the bearing and cryogenic flow system. The tasks 
described in this report do not represent the total analysis effort for the design 
and development of the BMT. The enclosed work was done to support the development of 
the BMT in specialized areas as problems occurred and to support the stated objec- 
tives. BMT test data are reduced, evaluated, and correlated with ana ,yses where 
applicable. Since test data applicable to the turbopump bearing system are limited 
and the system and operating conditions are of such _ complexity that theoretical 
modeling requires unprecedented extrapolations, additional data from the BMI are 
required to fully substantiate turbopump bearing analysis results and bearing life 
predictions. 
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2.0 SUMMARY 


Thermal and mechanical models of high-speed angular contact ball bearings 
operating in LOX and LN ? have been developed and verified with limited test data 
during the contract period. The work thus far demonstrates that a combined 
analytical-test approach is required to understand and evaluate these severe and 
complex bearing operating conditions. 

Significant progress has been made in understanding and adapting the SHABERTH 
bearing analysis computer program for evaluation of cryogenic bearing systems. This 
program was developed for the analysis of jet engine shaft/bearing systems operating 
above room temperature with normal hydrocarbon lubricants. Through continued 
analysis effort, it has been possible to adapt this tool to the evaluation of shaft 
bearing systems operating in cryogenics. 

Effects such as fluid drag, radial temperature gradients, outer race misalign- 
ments and clearance changes have been simulated and evaluated. In addition, the 
speed and preload effects on bearing radial stiffness has been evaluated. The 
SHABERTH program has also been used to provide contact stresses from which contact 
geometry has been calculated to support other analyses such as the determination of 
cryogenic fluid film thickness in the contacts and evaluation of surface and subsur- 
face stresses necessary for bearing failure evaluation. This program has been a 
vital tool for the thermal analysis of the bearing in that it provides the heat 
generation rates at the rolling element/race contacts for input into a thermal model 
of the bearing/shaft assembly. 

Although the SHABERTH computer code {See Reference 41) has a thermal subroutine, 
it is not capable of accommodating varying thermal properties and two phase flow. A 
bearing thermal model with this capability has been developed using the SINDA thermal 
analyzer. The bearing/shaft assembly and cryogenic coolant flow are thermally 
simulated by a nodal network and the model provides a detailed description of the 
bearing and assembly component temperatures and gradients. The model accounts for 
the various heat transfer regimes encountered in a two phase flow system and allows 
the assessment of parameters such as load, speed, flow rate, fluid subcooling, etc. 
on bearing component temperatures. Bearing contact surfaces are modeled in detail, 
to assess the maximum surface temperatures, and complete descriptions of component 
temperature gradients are provided. This can be important in assessing the charac- 
teristics such as fluid film thickness, film lubricant requirements, and effects on 
contact friction. Test data from the BMT has been used when possible to calibrate 
and verify the thermal model. Additional test data at varied tester operating 
conditions will be required for further model verification. 

Iteration between the SHABERTH bearing model and the SINDA thermal model is 
necessary for complete analysis of the bearing/shaft assembly. This is necessary due 
to the strong coupling between bearing component temperatures and internal operating 
clearances. As the inner race temperature exceeds that of the outer race, thermal 
growth reduces or eliminates the operating clearances. This in turn increases the 
contact stresses and forces which, increase the frictional heat generated. This is a 
cycle that will cause premature bearing failure if it proceeds unchecked. This 
condition has been simulated by utilizing the SHABERTH bearing program and the SINDA 
bearing thermal model. By this analysis method, preliminary bounds have been estab- 
lished for stable operation in LN 9 , These limits have been established in terms of 
fluid flow, fluid inlet temperatures, and axial load for a shaft speed of 30,000 RPM. 
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It is very important to conduct this analysis procedure to determine preliminary 
bounds for stable operation in LQX. Additional analyses of this nature have been 
done for a range of shaft speeds to determine the effects of shaft speed on bearing 
operating characteristics and temperatures. 

The importance of understanding and developing reliable analysis tools for esti- 
mating the mechanical/thermal interactions of the bearing system and these effects on 
bearing performance and life cannot be over emphasized. There are many deficiencies 
in the analytical methods available and applicable test data is nonexistent. There- 
fore, these gaps must be filled with additional BMT test and operational data before 
reliable predictions concerning Shuttle Main Engine turbopump bearing performance and 
life can be made. 

There is a strong mechanical /thermal interaction in the turbopump shaft/ bearing 
system. This is especially critical in cryogenieal ly cooled systems because of the 
drastically changing cooling ability as the fluid changes from liquid to vapor. As 
an example, an increase in load can cause transition from forced convection liquid 
cooling to forced convection vapor cooling with a significantly large increase in 
component temperatures . Such increases can rapidly eliminate internal clearances, 
increase contact stresses and frictional heat causing further heating until failure 
occurs. This unstable condition can further be aggravated by extremely high temper- 
atures in the contact areas which can cause breakdown of the dry film lubricant, 
significantly increasing the frictional heat generation. 

In addition to rapid bearing failure, less severe cases can cause failures 
different from . the subsurface fatigue failures experienced by conventional oil 
lubricated bearings. As will be discussed later, increased surface friction moves 
the maximum shear stresses to the surfaces and therefore, the failure occurs in the 
surface rather than below the surface. Furthermore, the high temperatures predicted 
will significantly degrade the material properties in the surface and this, combined 
with loss of lubricant, allows surface wear to become a potential failure mode. Due 
to the steep thermal gradients at the bearing contact surfaces, there is a strong 
possibility that thermally induced stresses are significant contributors to bearing 
failure. These conditions and failure modes are vastly different from those experi- 
enced by conventionally operated and lubricated rolling bearings. Therefore, the 
data base for conventional bearing life predictions has severely limited validity 
when extrapolated to high-speed bearings operating in cryogenics. 

Estimating coolant flow requirements for cryogenically cooled bearing systems, 
typical of the SSME turbopumps, requires consideration of mechanical /thermal coupling 
of the system. Analysis has shown that system heat generation is not a simple 
function of applied load. As previously discussed, the establishment of a tempera- 
ture gradient across the bearing will increase the frictional heat generated. For 
example, the coolant outlet temperature is essentially linear with flow for a con- 
stant heat input which, for a bearing of uniform temperature, would indicate a con- 
stant load. However, for a condition (this could be a low-flow case) that allows an 
adverse temperature gradient to establish itself across the bearing, the heat gen- 
erated will be considerably greater for the same axial load causing a higher coolant 
outlet temperature. Therefore, rather than a linear coolant outlet temperature in- 
crease with decreasing flow, the outlet coolant temperature increase can become 
exponential with decreasing coolant flow. The required coolant flow should be 
specified to maintain the essentially linear relationship between flow and outlet 
coolant temperature for given load and bearing speed conditions. 
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Contact stress analyses show that the maximum shear stress occurs in the surface 
of the contact for friction factors greater than 0.19. For marginally lubricated 
bearings operating in cryogenics, it is very likely that the contact friction factor 
is greater than 0.19. This is especially true when the bearing contact temperatures 
are significantly above the cryogenic saturation temperature. This condition pre- 
cludes any liquid film support between the contact. In addition, many dry film 
lubricants fail at temperatures approaching 500-700°F. Consequently, surface failure 
modes are strong candidates for bearings operating in these environments. 

Further research is needed to characterize the performance of dry film lubri- 
cants operating in cryogenics and the effects of contact surface temperature on the 
performance of these lubricants. The contact lubrication and traction, for bearing 
systems similar to those studied in this report, probably vary considerably with 
operation time and contact temperature. A key part of the bearing thermal analysis 
is the determination of the friction heat generated in the bearing contacts. This 
requires knowledge of local friction forces and slip velocities and the interdepen- 
dency of these variables for the specific lubrication system being analyzed. For 
example, the friction coefficient is generally a function of slip velocity and may be 
enhanced by metal oxide formation when operating in liquid oxygen. Additional test 
data characterizing these effects are required before reliable contact heat gen- 
eration estimates can be made. 

Additional work is also needed to characterize the heat transfer mechanisms for 
high-speed bearings operating in high-flow, subcooled, cryogenic systems. No data 
has been located in the literature search that is directly applicable to these 
conditions. 
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3.0 BEARING TESTER DATA EVALUATION AND TEST SUPPORT 
3. 1 LQX Test Incident Analysis 

In support of the Bearing and Material Tester (BMT) incident investigation, a 
scenario has been developed which describes the potential hazards associated with 
operating the tester bearings in the temperature regimes that produce two phase con- 
ditions in the coolant flow. The synergism of heating rate and flow pressure losses 
is much more pronounced in the two phase flow regime than in the single phase regime. 
A simplified flow network for the BMT is shown in Figure 3.1.1. Since each pair of 
bearings is cooled by flow circuits in parallel, the flow resistances of each circuit 
must be equal for equal flow and pressure loss. Furthermore, any flow losses up- 
stream of the bearings must be the same in each circuit. 

As indicated in the figure, this is not the case for the seal leakage. This is 
piobably not a serious factor since the absolute value of the real leakage is small. 
However, if the bearing component temperatures exceed the saturation temperature of 
the LOX coolant flow, vapor generation can significantly increase the resistance to 
flow through the hot bearing. This in turn can cause a reduction in flow, reduced 
cooling, and further increases in bearing temperature. 

The effect of increased vapor generation on flow pressure loss is shown in 
Figure 3.1.2. This is the correlation for steam/water mixtures developed by 
Martinelli and Nelson. Although the correlation is for steam/water flow, it is 
general in the sense that the same effect is experienced with other two phase mix- 
tures. As shown in the figure, there is a significant increase in pressure loss over 
the single phase value as the percentage of vapor is increased. In a parallel flow 
network similar to the BMT coolant circuit, this will reduce the flow to the bearing 
running the hottest and producing the most vapor. 

Another problem with operating in the two phase regime is the possibility of 
exceeding the critical heat flux for nucleate boiling. This can cause a rapid 
increase in component surface temperature. An example of this characteristic is 
shown in Figure 3.1.3. Shown in this figure is the measured outer race temperature 
for the failed bearing superimposed on the pool boiling curve for oxygen. The 
bearing temperature shown is for the outer race which is the lowest expected temper- 
ature for any of the bearing components (inner race, balls, etc.). The other bearing 
components will be running at a temperature in excess of the measured outer race 
temperature. 

Referring to Figure 3.1.3, it is evident that a small increase in load (heat 
flux) will cause the operating conditions to immediately shift to point E with a 
corresponding large increase in component temperature. The shift from C to D is not 
reversible. A reduction in load will cause the operating point to follow the path E 
to D and then move to the solid part of the curve horizontally across from point D. 
It is recognized that the boiling curve shown is for pool boiling and does not 
completely describe the high velocity flow through the bearings. In general, the 
velocity effect will be to increase the critical heat flux point. However, as stated 
previously, the measured temperature is the lowest expected for any bearing compo- 
nent. These two opposing factors will tend to cancel and the possibility of the 
operating point shifting as previously described is highly probable. 


Supporting data indicating possible effects of vapor generation in the load side 
coolant path are shown in Figures 3.1.4 and 3.1.5. This data represents the LOX 
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inlet manifold pressure for the flow circuits indicated. The mean steady pressure 
for the load side (P, 00 .) is about 440 psi, and the equivalent value for the drive 
side is about 419 psi. There is, however, a bias difference of about 7 psi at the no 
rotational point (time=0). Considering the bias, there is about a 14 psi difference 
between the two pressures with the load side being the greater. This could be caused 
by vapor generation in the load side circuit (bearings 1 and 2). 


Based on the above discussion, the following conclusions can be stated: 

(a) First order instabilities can occur when operating in a two phase regime 
near the peak flux region. This can cause a rapid increase in component temperature 
with a small increase in heat load. It appears, as a first order approximation, that 
bearing 2 was operating near this condition prior to the incident. 

(b) The parallel flow arrangement in the BMT requires stable flow resistances 
in both paths for equal coolant flow. Large variations in flow resistance can occur 
due to unequal heat generation in the bearings if the flow is two phase. 

(c) The combination of coolant flow reduction, increased vapor generation, and 
loss of cooling efficiency can cause a "run away" situation causing rapid increase in 
component temperature and subsequent failure. 


(d) Based on the above, the BMT bearing components should not be operated above 
the coolant, saturation temperature unless absolutely necessary. If operating above 
the saturation temperature is required, the operating condition should be evaluated 
with respect to the peak flux condition to assure adequate cooling margin and suffi- 
cient coolant flow to each flow circuit. 


3.2 Recommended Coolant Temperature Limits for 57 mm Bearing Operating in LNp 

Since the liquid density of LN ? is approximately 43% less than LOX at the nomi- 
nal operating temperatures and pressures for the BMT, the viscous work done on the 
LNp is estimated to be approximately 43% less than the work done on the LOX, provided 
the tester configuration was identical for the LN ? and LOX tests. The overall energy 
balances (Figure 4,3.6) for the LN,, and LOX tests show that the energy input to the 
coolant is about 11% higher for trffe LN„ tests. This apparent contradiction may be 
due to any of the following: L 

a) Changes made to the BMT following the LN, tests. Changes were made to the 
slinger and could account for part of tire apparent high heat input during 
the LN^ tests. 

b) The LOX may provide a better lubrication medium than the LN ? . Since LN ? is 
inert, any oxide films would not be reformed once worn avray. Such films 
could improve the friction coefficient. 

c) Inaccurate measurement of coolant temperature and flow. 

The conservative approach is assume that the tester changes will reduce the 
viscous work by the density ratio (.70) and therefore the heat generated in the bear- 
ing should be limited such that a maximum value of contact stresses (M25 K psi) will 
not be exceeded. The value selected corresponds to a 100°C temperature difference 
across the bearing, as shown in Figure 3.2.1. This is a relatively high value for 
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sustained. operation (For sustained running, the Hertz stress should be limited to 
^350 K psi) but should be adequate for establishing cutoff limits. 

Reducing the previously estimated viscous energy for LOX (Figure 4.3.6) to 
account for the lower UL density, gives a value for LN„ viscous work of 106.3 kW. 
The frictional heat generated by the bearings at 100°C AT is 15.2 kW. This gives a 
total heat input into the LNp of 121.5 kW. Based on the above values, the data 
provided in Figure 3.2.2 can ife generated. At flow rates less than 10 lbs/sec, the 
limiting criteria is to maintain the coolant outlet temperature about 30° below the 
saturation temperature. Since the outlet temperature for the second LN g tests (see 
Figure 3.2.3) rose about 30°F, restricting average inlet flow to 160°R vfould give an 
outlet temperature of 190°R which is approximately 30° below the saturation tempera- 
ture of LN« at 400 psi a. This provides a safety factor of 2 with regard to the 
outlet flovr conditions. 


The above values provide limits based 
necessary to monitor and limit the energy 
lish the limiting AT across the bearing, 
were estimated assuming a 100°C radial 
(limits Hertz stress to M50 K psi). The 
bearing pair was estimated to be 32 kW. 

Figure 3.2.4 were generated. The curves 
zontally (dashed lines) at specific values of 

available AT is approximately 60°R (T_ at - T rw ,. ueyrt^ ut mis vaiue was 

used as a safety factor. At a flow rate of lOlbs/sec (total), about 12 degrees AT 
occurs across the bearing. This would allow about 18° for viscous work on the fluid 
excluding that done by the bearing. This is the minimum flow that will meet the 
imposed temperature restrictions. However, since a new tester configuration is being 
tested for the first time, flow rates less than those previously used, 12.6 Ibs/sec, 
are not recommended. Based on a flow of 12.6 Ibs/sec, an inlet temperature of 160°R, 
a ^'’ting AT of 18°R across the tester, the cutoff outlet temperature should be 
1/8 R or -282°F, and the AT across the bearing should not exceed 9°R. 


on an overall energy balance. It is also 
generated by each bearing pair. To estab- 
the bearing viscous work and friction heat 
temperature difference across the bearing 
combined viscous and friction heat for one 
Based on this value, the data provided in 
in Figures 3.2.2 and 3.2.4 are cut hori- 
aT available across the tester. The 
). Thirty degrees of this value was 


To gain a clearer understanding of the magnitudes of the viscous and bearing 
frictional heat generation and to minimize the risk to the BMT, it is recommended 
tha + a low speed (^10,000 RPM) test be performed with no applied load prior to runs 
at^design speed. Data from this test could be used to verify heat generation esti- 
mates for updating temperature cutoff values. In view of the uncertainties in the 
analysis, the cutoff values provided should be considered preliminary for design 
speed and load tests. These analyses should be verified and/or updated based on 
additional data from the recommended low speed test. These cutoff values should be 
adequate for the low speed tests. In addition, it is recommended that the outer race 
temperatures be limited to 200°R. This is about 20°R below the saturation tempera- 
ture at 400 psia. 


3 -3 Comparison of Bearing Tester Data with Thermal Model Results 

A comparative analysis was performed between the test data for an LI\L run and 
temperature data from the thermal model for various operating conditions. “The outer 
edge of the outer race is divided into several nodes as shown in Figures 3,3.1 thru 
3.3.4. Model temperatures are shown superimposed on the outer race thermocouple 
measurement. The operating conditions for the thermal analysis are noted on each 
figure. 
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Figure 3,3.1 shows the data from thermocouple number T1004 located on the back 
side of the outer race of bearing number 1. The dashed lines indicate the predicted 
thermal model temperatures for the nodes illustrated in the sketch,. As noted in 
Figure 3.3.1, this analysis was conducted for an LN ? coolant flow rate of 6.3 lb/sec 
per bearing set, an inlet coolant temperature of -3x)5 0 F and an axial reaction of 500 
pounds. 


These analyses considered an axial load on the carrier of 2,500 lbs caused by 
the pressure differential across the carrier. This loads bearings 2 and 3 to 2,500 

«A„ a , n . d x unloads bearings 1 and 4 to 500 lbs (nominal preload for all bearinqs is 
1,000 lbs). a 


A comparison of these data does not indicate a significant flow unbalance 
between bearing pairs. There does appear to be measurement problems with T1004; the 
outer race temperature of Bearing 1. This temperature rises only about 4° during the 
test, while the comparable temperature of bearing 4 rises 17°F. Further comparison 
of bearing 1. outer race temperature rise (T1004) with comparable measurements taken 
during rotational test 03-30-84, which reached a speed of about 15,000 RPM indicates 
that the measurement is suspect. The outer race temperature rise for bearing 1 
during the 03-30-84 test rose 3.3°F which is almost as much as the data shown in 
3.3.1. Since the speed during this run was approximately 30,000 RPM the outer race 
temperature rise should be significantly higher. Assuming T1004 is a bad measure- 
ment, the thermal model results compare reasonably well with the measured data for a 

o n L fl ° W u ° f 6-3 Tbs/ ? e ? P er bearing pair. These results are shown in Figures 
3.3.2 through 3.3.4, Additional tester data is obviously necessary to arrive at firm 
conclusions concerning bearing loads coolant flow effects. 

3 -4 Analysis of Bearing Temperatures and Operating Characteristic s 

3 *^*^ Sensit ivity of Node Temperatures as a Function of Flow Rate and Degree of Sub- 
cool ing 


. , da ^ a ,base developed as a part of the bearing thermal model was set up 
with heat generation conditions to simulate a 500 pound axial reaction on bearing 1 
and a 2500 and 6000 pound axial reaction on bearing 2 for a 57 mm bearing set. 
During the course of this investigation four different flow rates thru the bearing 
set; 3 *2, 6.3, 9.6 and 12.8 Ibs/sec were considered. For each different mass flow 
rate and loading condition, the thermal model was exercised for two different values 
of subcooled temperature; -285°F and -305°F. The frictional heat generated was 
determined based on a friction coefficient of 0.2 and a AT across the bearinq as 
predicted by the SHABERTH bearing analysis program. 


The sensitivity of several node temperatures was examined to determine the 
effect of changing mass flow rate, degree of subcooling and loading conditions. In 
figure 3.4.1 a plot of the temperature of the hottest track node on the inner race 
the mass flow rate of the coolant for several loading conditions is shown for 
different degrees of subcooling. 


As can be seen in Figure 3.4.1 for bearing 1, the effect of coolant flow, and 
suDCooling on race temperature is about the same for the range of variables inves- 
tigated. The higher load cases, 2500 and 6000 pounds on bearing 2, shows a stronger 
influence on coolant flow on track temperatures; especially for the case of coolant 
inlet temperatures of -305°F. The bearing 2 data in Figure 3.4.1 show that for the 
lower flows, the race temperature becomes less sensitive to subcooling. It also 
shows that for marginal subcooling, increasing coolant flow is not as effective in 
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control! ins r ace temperature. Consequently, both flow and subcooling need to be 
adequate for proper bearing cooling. 

The dependence of bearing cooling on coolant flow is complicated due to the 
dependence of bearing internal clearances on component temperature. As internal 
clearances are reduced due to thermal growth of the inner race, with respect to the 
outer race, the normal forces and stresses at the contact surfaces are increased 
resulting in further increases in frictional heat generation and increased tempera- 
tures. Therefore, the bearing is not necessarily a constant heat generator for fixed 
load, speed, and varying coolant flow. Local subcooled boiling is accounted for by 
using two-phase boiling heat transfer correlations. This represents the physical 
system provided the outlet coolant temperature does not equal or exceed the coolant 
saturation temperature. As shown in Figure 3.4.2, for a pressure of 400 psia, the 
saturation temperature of LN, is -240°F. For a inlet coolant temperature of -285°F 
and an axial reaction of 250Cr and 6000 pounds, the bearing outlet coolant temperature 
approaches the LIn 2 saturation temperature at low flows. Outlet pressure slightly 
less than 400 psia could easily result in a condition of bulk vapor generation at the 
bearing outlet, increased flow resistance, and subsequent bearing failure due to 
inadequate cooling. 

Figures 3.4.2 and 3.4.3 illustrate the dependence of frictional heat load, 
bearing outlet temperature, and coolant flow for the fixed conditions of axial load 
and coolant inlet temperatures. The 6000 lb axially loaded case shows an increase in 
frictional heat load as the coolant flow rate is reduced. As the coolant flow is 
decreased, there are two factors causing an increase in bearing outlet coolant 
temperature: increased heat load and flow reduction. If thermal equilibrium cannot 
be attained at a higher temperature compatible with bearing survival then failure 
occurs. The 2500 lb load case did not experience a heat load increase over the flow 
range investigated for either -285°F or -305°F inlet coolant temperature (F : gures 
3.4.4 and 3.4.5 respectively). 

Test data from the planned BMT tests are critical to guide and verify the 
analysis effort. This program is providing considerable insight into the thermo- 
mechanical interactions of the bearings which will, when verified, provide extremely 
valuable information for improving the design and life of flight hardware. 

3 - 4 - 2 Bearing Operating Characteristics as a Function of Speed and Axial Load 

Observed wear tracks of bearing races and rolling elements can be used as a tool 
for estimating the loads experienced by the bearing components by matching the 
measured contact angles and wear tracks with predicted values. Shown in Figures 
3.4.4 through 3.4.7 are contact angles and wear track half widths (semi -major axis of 
contact ellipse) as a function of shaft speed and axial loads. Future work in this 
area should include the correlation of test observations from the BMT with the 
enclosed data generated from the SHABERTH bearing code. 

3 - 4 ’ 3 Bearing Coolant Flow Effects on Bearing Heat Generation, Component 

Temperatures, Operating Clearances, Hertz Stresses, and Contact Angles 

The 57mm angular contact ball bearing (L0X turbopump turbine end bearing) was 
analyzed to investigate the effects of varying the coolant flow rate (LN„) on bearing 
temperatures and operating characteristics such as clearances, contact angles and 
contact stresses. These effects were investigated for an LN„ coolant inlet tempera- 
ture range of -285° F to -305° F. 2 
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The bearing thermal model provides a detailed representation of the bearinq 
component temperatures and gradients. Due to high frictional heat generated in the 

u 0f * he b 1 earin 9 and the 1ow (cryogenic) temperatures of the coolant, 
y 4 .u hl9 ! ) i tber L mal gradients are present in the contact areas. Figure 3.4,8 
provides the total heat generated by a bearing pair and the coolant outlet tempera- 
ture as a function of the coolant flow rate for the conditions indicated. This 
pe r forme : d t°. r an axial reaction on bearing #2 of 2500 lbs. The heat 
? “* as | h ?wn in Figure 3.4.8, includes the viscous energy dissipated to the 
fluid and the friction heat generated in the bearing contacts. As indicated, there 
is very little change in slope of the heat generation curve as the coolant flow 
decreases. This implies that the net heat, due to the loading conditions indicated, 
IL, \ large enough to significantly reduce clearances and contact angles as the 

fl ° W de£ T eases - Therefore, thermal runaway should not be a problem for these 
operating parameters. r 

. ^ S °' indicated in Figure 3.4.8, the coolant outlet temperature is approach- 

ing the coolant saturation temperature. If the coolant leaving the bearing pair 
becomes saturated, . net vapor generation occurs and vapor effects on coolant pressure 
losses can become important considerations. P 

•s n Q F a! 1 H tl ? n ? 1 in e | t r tl e profiles over the contact surfaces are shown in Figures 
r\ . and 3.4.10 for the inner and outer races, respectively. The detail of the 
bearing thermal models used for the analysis allows this estimation of the distribu- 

d ?trih.iHnn C c ^ W \ th l n ^ el11ptical contacts and the resulting temperature 
? ver . the contact ellipse area. The frictional heat rates obtained from 

sJrf/r H p AB n^ T « ^ ode con L tacts is distributed to the corresponding 

surface nodes in the thermal model to obtain temperature profiles for the bearinq 

rnnf°^ ntS f Figure 3 * 4 / 11 shows the maximum temperatures estimated in the bearing 

temner^iirtr ^ h ^ on h dltl ?P s sb ° wn - As expected, the estimated maximum contact 
temperatures are higher than the coolant saturation temperature, therefore, the 

Sc a - 1 H?rt5I )nt i CtS are ? pera J 1n 9 in 9 as rather than liquid. Consequently, the LN 9 acts 
as a medium for removing the heat generated and little or no lubrication is provided. 

ft..**/™™?®™ 9 ® beari ( , 9 component temperatures (rolling element, inner race, and 
sharfpth t at ^ e show " ™ Figure 3-4.12. This information is necessary input to the 
S??::, bearin 9 model and is obtained from a code incorporated into the SINDA 
thermal analyzer program which averages individual node temperatures taking into 

of over eon nndlT ° f - + aCh T* 6 * A * th . 1s time > the bearing thermal model is composed 

tures indiV^tP^ n A S ’ K Vi 9U k re 3 /. 4 ’! 2 .* u which shows the average component tempera- 
onfpr’ th - at the bal ] has the highest average temperature with the inner and 

race following. Considerable vapor will be generated locally due to the 
component temperatures that exceed the LN ? saturation temperature. Vapor that is 

nutfp^tpml 009 ! y W1 J CGndensed and 1 ^ ave the tester as liquid since the bulk 
dieted Vn ^ngure^ °4 8^* C °° 1ant remains be1ow the saturation temperature as in- 

chn,m Th i e n l°- SS ° f or °P eratin 9 clearance as a function of coolant flow is 

2SJ" in 3 ' 4 ; 13 a1 f9 «1th increase in contact stresses as the flow is 

hPi^w ahniS? n h0 i W K/ the reduGtl °! 1 in clearance becomes more pronounced at flow rates 
\ 8 1b (t ec pe , r bearin 9 P air * Figure 3.4.14 illustrates the trend in 
contact angles as the coolant flow is decreased. As indicated, the change in contact 

250C> e aS d fUnCtl0n 0f the CGolant flow rate is very small for an axial reaction of 
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A similar analysis was conducted for an axial reaction of 6000 pounds on bearing 
#2. the results of this analysis are found on Figures 3.4.15 through 3.4.19. Figure 
3.4.15 provides the total heat generated by a bearing pair and the coolant outlet 
temperature as a function of coolant flow for the conditions noted. The heat gen- 
erated includes the viscous energy dissipated to the fluid and the friction heat 
generated in the bearing contacts. Test data from the BMT have been used as an aid 
in estimating the viscous energy dissipated in the LN 2 coolant. Additional data are 
necessary to improve the confidence in these estimates, especially for LOX cooled 
bearings. As shown, there is a change in slope of the heat generation curve as the 
coolant flow decreases below about 8 Ibs/sec per bearing pair. This is the result of 
the "amplification" of loads due to reduction in clearances and contact angles as the 
coolant flow is decreased. It is also significant to note that the coolant outlet 
temperature is approaching the coolant saturation temperature. In the event the 
coolant outlet becomes saturated, net vapor generation occurs and vapor effects on 
bearing coolant pressure losses can become significant. 

As was previously done, frictional heat rate distributions from the SHABERTH 
bearing model were input into the SINDA thermal model to provide temperature 
distributions. Frictional heat rate profiles over the contact surfaces are shown in 
Figure 3.4.16 and 3.4.17 for the inner and outer races. These heat rates are 
distributed to the proper elements in the thermal model to obtain temperature 
profiles for the bearing components. Shown in Figure 3.4.18 are the maximum 
temperatures estimated in the bearing contacts for the conditions noted. As indicat- 
ed, the estimated contact temperatures are relatively high compared to the cryogenic 
coolant temperature. Since the saturation temperature of LN 2 is ~240°F at 400 psia, 
the bearing contacts are operating in gas rather than liquia, and liquid film sepa- 
ration of the contacts does not exist as in conventionally lubricated bearings. 
Consequently, the LN„ is expected to provide very little if any contribution to 
bearing lubrication and serves only as a medium for removing the heat generated, as 
previously discussed. This emphasizes the importance of the dry film transfer 
lubrication technique employed in this bearing. Loss of lubricant can cause in- 
creased friction in the contacts increasing not only wear but heat generation which, 
due to increased temperature, can cause loss of operating clearances and ultimately 
bearing failure. 

As discussed, the SINDA bearing thermal model is a finite element representation 
of the bearing components and provides a two dimensional representation of the 
temperature profile in the bearing. The SHABERTH model, however, requires the 
average temperature of each rolling element component (inner race, ball and outer 
race} to be supplied as input. To meet this need, a temperature averaging capability 
for each element of the bearing has been developed in the thermal model. These 
average component temperatures are shown in Figure 3.4.19. The ball has the highest 
average temperature with the inner and outer race following. The ball and inner race 
average temperatures exceed the LN 2 saturation temperature for the conditions noted 
and, consequently, considerable vrpor will be generated locally. Since the bulk 
outlet temperature of the coolant remains below the saturation temperature, this 
vapor will be condensed and leave the tester as liquid. It is also important to note 
that the inner race and ball temperature exceed the outer race temperature. This 
causes a loss in bearing operating clearances due to the thermal growth of the ball 
and inner race relative to the outer race. 

The loss of internal or operating clearance as a function of coolant flow is 
shown in Figure 3.4.20 along with the increase in contact stresses as the coolant 
flow is reduced. As indicated, the reduction in clearance becomes more pronounced at 
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flow rates below about 8 Ibs/sec per bearing pair. Figure 3.4.21 shows the reduction 
in contact angles as the coolant flow is decreased. This is the principal contribu- 
tion to the increase in contact stresses. These analyses indicate that bearing 
failure can occur if conditions are present that allow thermal growth to reduce the 
operating clearances such that additional heat is generated, thus further reducing 
the clearances. 

3.4.4 Bearing Viscous and Frictional Heat Generated as a Function of Speed and 
Radial Load 


Parametric analyses were done to determine the relationship between the heat 
generated, in a bearing pair, and shaft speed and radial load. The SSME LOX pump 
turbine end bearing pair was modeled on the SHABERTH computer code. Heat is generat- 
ed during bearing operation by two mechanisms; viscous energy dissipation into the 
coolant, and friction generated in the bearing contacts. The SHABERTH code is not 
used to estimate the viscous energy. Although SHABERTH provides a capability to 
estimate the viscous energy dissipated due to ball drag, this is not a complete 
description of the energy dissipated to the fluid. For example, fluid frictional 
energy is dissipated, due to the relative velocity of the fluid and shaft, as the 
coolant flows from the first to the second bearing. In addition, the rolling ele- 
ments not only orbit the shaft but spin relative to their own axis. This spin motion 
produces a surface velocity relative to the fluid velocity contributing to the 
dissipation of viscous energy. 

Figure 3.4.22 is a representation of the viscous energy dissipated as a function 
of shaft speed for a pair of bearings. Shown in Figure 3.4.23 is the friction heat 
generated in the bearing contacts as a function of radial reaction and shaft speed. 
Figures 3.4.24 and 3.4.25 show the frictional heat generated at the inner and outer 
races as a function of shaft speed and radial load. Figure 3.4.26 is a cross plot of 
the data shown in Figure 3.4.23. It is evident that the heat generated is nonlinear 
with shaft speed. Since the LOX turbopump coolant flow increase, with shaft speed, 
is essentially linear, an increase from rated power level (RPL) reduces bearing 
cooling margin as power level increases. 

Similar parametric analyses were performed to determine the relationship between 
the friction heat generated in a bearing pair and shaft speed and axial load. Figure 
3.4.27 is a representation of the friction heat generated in the bearing contacts as 
a function of axial reaction and shaft speed. Figures 3.4.28 and 3.4.29 illustrate 
the friction at the inner and outer races as a function of the shaft speed and axial 
load for the conditions noted. Figure 3.4.30 is a cross plot of the data shown in 
Figure 3.4.27. It is apparent that the heat generated is nonlinear with shaft speed. 
The friction heat generation data presented assumes no thermal gradients across the 
bearing and, therefore, should be considered as trends only. A combined ther- 
mal/mechanical analysis will be required to accurately estimate the friction heat 
generated especially at the higher loading conditions. 

3.4.5 Investigation of Bearing Temperature Gradients on Operational Characteristics 

A single 57 mm ball bearing was modeled on the SHABERTH computer program. 
Several temperature gradients were imposed and their effects on contact angle, normal 
forces, contact stresses, operating clearances, and heat generation were determined. 
The results, shown in Figures 3.4.31 through 3.4.34 are for a bearing with .0036 inch 
radial play at a contact angle of 20°. 
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Figure 3.4.31 Shows the variation of inner and outer race contact angles and 
stresses as a function of bearing temperature difference. The positive temperature 
gradient is from the shaft to the bearing housing. The shaft and inner race are 
assumed to be at the high temperature and the ball temperature is assumed to be the 
average of the inner and outer race. As shown, the operating contact angles exhibit 
a significant decrease at the higher temperature differences with a corresponding in- 
crease in contact stresses. Figure 3.4.32 shows the effect of the temperature 
differences on the bearing operating clearances. As the clearances go from positive 
to negative, an internal preload is introduced which in effect increases the bearing 
preload. As shown, this occurs for a rather small temperature difference. It should 
be pointed out, however, that these analyses were made using the worst case stack up 
of tolerances to minimize the bearing clearances. 

Shown in Figure 3.4.33 are the normal forces per ball for the inner and outer 
races as they are influenced by the bearing temperature differences^ The forces 
appear to be approximately linear up to a temperature difference of 150 C and become 
a power function or exponential function as the AT increases. Shown in Figure 3.4.34 
is this effect on the increase in heat generated at the contact surfaces. After 
about 200 I> C aT, there is a sharp rise in heat generated as the AT across the bearing 
is increased. The increase in heat generated will increase the bearing AT and, in 
turn, the AT effects will increase the heat generated. This can result in an unsta- 
ble condition causing premature bearing failure. 

This analysis was repeated for a bearing having zero radial play at a contact 
angle of 20°, which is representative of the test bearing. These results are shown 
in Figures 3.4.35 through 3.4.38. In comparing Figure 3.4.35 with Figure 3.4.31, it 
is evident that eliminating the radial play reduces the contact angle and increases 
the contact stresses, as should be expected. Comparison of the operating clearances 
(Figure 3.4.36) indicates a slight reduction in clearance for the case of no radial 
play. This results in an increase in frictional heat generated at the contact 
surfaces as can be seen by comparing Figure 3.4.38 with Figure 3.4.34. 

Elimination of the radial play slightly increased the sensitivity of operating 
conditions to bearing temperature differences. The previously discussed charts are 
labeled "minimum clearances". However, there may be other combinations of clearances 
that produce slightly tighter operati no clearances. The lower range of operating 
temperature for these analyses was -18(rC or 168°R, which is representative of. the 
LH„ test environment. The higher temperature of the LOX (193 R) will have a slight 
effect on the clearance values. 

3.5 Pretest Predictions of Thermal Conditions for the BMT 


Estimates have been made for coolant and bearing temperatures expected during 
the first rotational test of the BMT (Build #1). The speed-flow profile for the 
first test is shown in Figure 3.5.1. The various conditions analyzed are shown in 
Figure 3.5.2. A total of 114 SHABERTH and SINDA computer runs were necessary to 
evaluate the conditions shown. 

The rate of energy dissipation due to viscous mixing was estimated for each 
value of shaft speed. Bearing friction heat was estimated for each load and shaft 
speed. This information is shown in Figure 3.5.3. These data were used in the 
bearing thermal model to estimate the thermal conditions of the coolant and bearing 
temperatures at the measurement locations in the BMT. 
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Temperatures for a coolant inlet temperature of -30Q°F and axial reaction of 
2,500 lbs are shown in Figure 3.5.4. A shaft speed of 25,000 RPM obviously produces 
the highest temperatures. The bearing outlet coolant temperature for this condition 
is -282°F. If a 30°F margin is added, the allowable saturated outlet temperature is 
-252°F and the corresponding saturation pressure is 290 psia. Therefore, the pres- 
sure at the bearing exit should be greater than 290 psia to meet the stated outlet 
temperature criteria. This will allow a pressure loss of 110 psi from tester inlet 
to bearing outlet assuming an inlet pressure of 400 psia. This condition is shown in 
Figure 3.5.5. As indicated, conditions at the lower speeds are less severe and the 
allowable pressure loss though the tester increases as speed decreases. Similar data 
for an axial reaction of 4,000 lbs are shown in Figure 3.5.6. As indicated in Figure 
3.5.5, the 4,000 lbs load and 25,000 RPM case allows a pressure { inlet-to-bearing 
exit) of about 100 psi. 

Two load conditions were evaluated for an inlet coolant temperature of -280°F. 
Results of the 2,500 lb case are shown in Figure 3.5.7 and Figure 3.5.8 provides 
similar results for a 4,000 lb axial reaction. As shown at 25,000 RPM, the coolant 
temperature at the bearing outlet is -262°F. Adding 30° of subcooling gives an 
allowable outlet of -232°F. Since the critical temperature of LN ? is -232.8 at 493 
psia, the above allowable outlet temperature exceeds the critical temperature. This 
requires the inlet and outlet coolant pressure conditions to be super critical. 
Since this is undesirable, an inlet temperature of at least -300°F should be provid- 
ed. The inlet pressure requirement is a function of pressure loss from tester inlet 
to bearing exit. A preliminary estimate of the pressure loss from inlet to the 
bearing exit is on the order of 100 psi. Since pressure loss estimates for the BMT 
coolant contain unknown error margins, the bearing inlet pressure should be as high 
as practical and the pressure losses during the low speed runs should be evaluated to 
determine the possibility of exceeding the exit temperature criteria at the higher 
speed. 

The outer race and ball track temperatures for bearing No. 2 are shown in Figure 
3.5.9 for loading and coolant inlet conditions indicated. This information was 
generated from the SHABERTH bearing code and the SINDA thermal model. As indicated, 
these temperatures are not expected to be severe for the load and speed conditions of 
the test. Shown in Figure 3.5.10 are the estimated average temperatures for the 
ball, inner race and outer race. 
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FIGURE 3.5.9 BALL TRACK AND BACK SURFACE OUTER RACE TEMPERATURES VS SHAFT SPEED 


OPERATING CONDITIONS 

COOLANT INLET TEMP - 300°F 
COOLANT - FLOH 6 Lb/8GR PAIR 
AXIAL REACTION (Lbs) 

— 4000 


•" BALL TRACK 




-100 ) — - — I r- 

IS IT 19 



— 1 OUTER RACE 

I BACK 

1 " 'I 

23 2S 


SHAFT SPEED (RPM) XIO 3 


FIGURE 3.5.10 


AVERAGE BEARING COMPONENT TEMPERATURE VS SHAFT SPEED 


OPERATING CONDITIONS T~ 

COOLANT INLET -300°F j 

COOLANT FLOH 6 Lbs/sec BGR PAIR 
AXIAL REACTIOH (Lbs) *- J 

— 4000 
2500 



32 


COOLANT 4T 


spectra research systems 


4.0 BEARING THERMAL MODEL DEVELOPMENT 


The development, refinement, modification and preliminary verification of a 
thermal model using geometrical configurations of the 57 mm LOX Bearing and Materials 
Tester has been accomplished during the contract period of July 1982 to April 1984. 
The primary analytical tools used in this development were: 1) SINDA - a numerical 
differencing computer code used primarily for thermal systems, 2) SHABERTH - a 
state-of-the-art bearing analysis program. 

4.1 Bearing and Materials Tester Data Base 

A data base composed of a finite element representation of a pair of 57 mm 
back-to-back duplex mounted ball bearings as well as the shaft, carrier, housing, and 
coolant has been established using the dimensions of the Bearing and Materials Tester 
(See Figure 4.1,1). In lieu of modeling the entire bearing, a pie-shaped slice as 
shown in 4.1.2, encompassing 1 of the 13 rolling elements in the bearing, was used as 
a boundary for the thermal model . 

Figure 4.1.3 illustrates the detail involved in the node analysis of the races 
and rolling element. This level of detail is necessary, especially in the track 
area, due to the large thermal gradients, present as shown in Figure 4.1.4 for the 
rolling element. This data indicates a rather large temperature gradient of 2300°F/- 
inch in the rolling element. 

A data base has been established for both LOX and LN^ coolant flow through the 
Bearing and Materials Tester (See Appendix A). This data base includes all node and 
conductor data along with temperature varying thermal property data and film coeffi- 
cient data. The data base has been designed to facilitate modification due to 
changing operating parameters such as flow rates, axial load, and degree of subcool- 
ing or inlet coolant temperature. The bearing thermal model is composed of 630 
diffusion nodes, 9 arithmetic nodes and 3 boundary nodes for a total of 642 nodes, as 
well as 1387 linear conductors. 

SINDA, a software system which possesses the capabilities that make it well 
suited for solving lumped parameter representations of physical problems governed by 
diffusion-type equations, is used in conjunction with the data base. The system is 
designed as a general thermal analyzer accepting resistor-capacitor network represen- 
tations of thermal systems. 

4.2 Bearing Thermal Model Temperatures 

The bearing thermal model, when executed, achieves a system heat balance for 
steady state conditions and when the model meets convergence criteria concerning the 
allowable temperature change between iterations, the final node temperatures are 
printed as shown in a representative case in Appendix B. Representative bearing 
component temperature distributions are shown in Figures 4.2.1 through 4.2.8 for 
conditions of loading, and coolant flow rate and temperature as noted. As shown, the 
gradients in the track area are significant for the higher loading conditions. Track 
surface temperatures, for the larger loads, are of concern due to the possibility of 
high temperature breakdown of the solid film lubricant. Since the track temperatures 
are critical to lubricant performance and in order to develop a better understanding 
of the effects of a distributed heat load in the contact area of the bearing, modifi- 
cations were implemented into SHABERTH in order to determine a node-by-node heat 
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FIGURE 4.1.3 DETAILED NODAL REPRESENTATION 


OF 57mm BEARING 



36 



spectra research systems 


ORIGINAL PAQF r- 
OF POOR quaCit? 


FIGURE 4.1.4. TEMPERATURE -VS- DEPTH BELOW SURFACE OF ROLLING ELEMENT 
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FIGURE 4.2.1 

ROLLING ELEMENT TEMPERATURE DISTRIBUTION FOR LOX PUMP TURBINE END BEARING 
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FIGURE 4.2.4 INNER RACE TEMPERATURE DISTRIBUTION FOR LOX PUMP TURBINE END BEARING 
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FIGURE 4.2.5 

ROLLING ELEMENT TEMPERATURE DISTRIBUTION FOR LOX PUMP TURBING END BEARING 





42 









o o 

-n X 





spectra research systems 









FIGURE 4.2.7 OUTER RACE TEMPERATURE DISTRIBUTION FOR LOX PUMP TURBINE END BEARING 
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generation value as a function of loading for input into the bearing thermal model. 
The modifications to 5HABERTH are further discussed in Section 4.5. 

Experimental rotational data from the BMT (Bearing and Materials Tester) is 
essential in order to verify the accuracy of the bearing model which has been devel- 
oped. 


4.3 Thermodynamics of the Bearing System 
4.3.1 Energy and Flow Balance Estimates 

Energy and flow balance analyses were conducted to investigate the differences 
in energy flow into and out of the tester during the LN ? tests and LOX tests. In 
addition, inlet and outlet pressure measurements were evaluated in an attempt to 
compare tester flows during the LN ? and LOX tests. Four test conditions were evalu- 
ated: (1) LOX flow, (2) LOX rotational, (3) LN 2 flow, and (4) LN ? rotational. The 
flow tests were evaluated to judge the validity^of pressure and temperature readings 
and to estimate the heat gain of the tester flow from the atmosphere. Figures 5.3.1, 
5.3.2, 5.3.3, and 5.3.4 provide a summary of the four test conditions. These figures 
show the test values and schematically locate the instrumentation of interest. Test 
identification information is shown in the upper right corner of each figure. 

The flow test data (tester not rotating) are shown in Figures 5.3.1 and 5.3.2. 
As indicated in the lower left corner of Figure 5,3.1, temperature measurements T.™, 
through T. Q03 are resistance bulb (RB) thermometers and the other temperature mea- 
surements ^re thermocouples. There appears to be a problem with the RB's during the 
LOX test since they show a decrease in temperature across the tester in the flow 
direction (see Figure 5.3.1). They also read about 6°F lower than the thermocouples 
in comparable locations. The RB's are in fair agreement with the thermocouples 
during the LN, flow tests (Figure 5.3.2). Since the RB's appeared to be in error 
during the LOK tests, the thermocouple measurements at the entrance and exit mani- 
folds were used to estimate the heat flow. 

Pressure measurement Pg is inconsistent with the upstream and downstream pres- 
sure measurements for both lOX test conditions. T, nn2 is not provided for the LOX 
tests and T^ is a bad measurement. As stated above , the LOX flow test data were 
evaluated to determine the ambient heat flow into the tester. As shown in Figure 
5.3.1, there is very little temperature gain as the flow moves from the inlet to the 
exit manifold. Using the stated LOX flow rate of 9.1 Ibs/sec. and the inlet and 
outlet temperature difference, an ambient heat input of about .5 Btu/sec was cal- 
culated from the LOX flow test data. 

The BMT thermal system was defined as follows: 
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where: 

m 

mass flow rate 


T 

temperature 


P 

pressure 


Q 

heat rate 


w 

power 


h 

fluid enthalpy 

subscripts: 

i 

inlet conditions 


0 

outlet conditions 


An energy balance on the system defined above is as follows: 

w Q = m(h 0 - h.) + || 
dE 

where: is the time rate of change of internal energy assumed to be zero, thus 

implying steady state conditions. 

The power term (to) is determined from the torque and speed measurements and, as 
stated previously, the ambient heat flow term was evaluated from the flow test data 
and assumed to be the same for the rotating tests. Shown on Figure 4.3.5 are the 
terms of the above heat balance for the four test conditions. As shown in the fig- 
ure, there is an excess of energy not accounted for in the steady state energy 
balance. The LOX test gives about 21 kw and the LhL test about 9.4 kw excess energy 
input into the system. This is a maximum error of 12% and 5% for the LOX and LN 2 
heat, balances. L 

Although this is fairly good agreement, there are several possible explanations 
for the energy imbalance. The calculation of power input to the tester is question- 
able because the torque is measured at the gear box input. This requires the gear 
box torque to be determined at no load conditions and subtracted from the test value. 
Since the loaded gear box torque is very likely greater, and the gear box torque can 
be affected by gear box temperature, etc., this method of determining the BMT torque 
has several uncertainties. A torque measurement at the gear box outlet is needed to 
eliminate these questions. 

Other sources of potential error are the temperature and flow measurements, and 
the storage of energy in the tester hardware. Since the tester housing temperature 
does increase during the rotation tests, some energy is being stored. An accurate 
estimate of the amount is difficult due to the complicated geometry of the tester. 
Assuming for the moment that the energy storage term is small and the power input is 
correct, it is interesting to note that the energy input to the coolant fl'ow (mAh) is 
about 11% greater during the LNp runs as compared to the LOX tests. Because of the 
higher LOX density it would be expected that the viscous work done on the LOX would 
be considerably greater than done on the LN ? . This could suggest higher bearing 
friction in the L^. L 

Based on the preceding evaluation, the following conclusions can be stated: 

1) An overall system energy balance is a valuable tool for evaluating the 
validity of the temperature, pressure, torque, speed, and flow data. 

2) To improve the capability of this technique several changes should be made 
in the tester instrumentation as suggested below: 
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FIGURE 4.3.5 ENERGY BALANCE SUMMARY 
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a) Provide a torque measurement at the drive end of the BMT. 

b) Add additional temperature and pressure measurements at the entrance 
and exit manifolds of the tester to assure more than one measurement 
at these critical locations. The temperature measurements should be 
that of the fluid and not the metal. 

3) It is highly desirable to determine the flow rate to each pair of bearings 
as well as the total flow through the tester. Since flow is such a crit- 
ical parameter, redundant measurements need to be considered. 

4 * 3 * 2 ° f Tester Int ernal Power Dissipation with Bearing Temperature 


To determine the sources and magnitudes of heat generated in the BMT, each 
component was evaluated and an estimate of the amount of energy dissipated was made, 
me total amount dissipated was then compared to the total energy transferred to the 
LOX coolant flowing through the bearings. As will be shown, the total enerqy es- 
timated from the component analysis is less than the energy transferred to the LOX 
coolant. The LOX coolant and component energy can be balanced, however, by assuminq 
a temperature gradient across the bearing. The temperature gradient causes increased 
internal loading and heat generation. 


The power estimates for each component in the BMT are shown in Figure 4.3.6. 
This information was calculated using semiempirical relationships provided in Refer- 
ence 2 and viscous ball drag obtained from the bearing computer program. Several 
assumptions were necessary to complete the power estimates. The most significant are 
as follows: 


1} The fluid has no rotational speed prior to entering bearings 1 and 4. 

2) Front faces (face to shaft shoulder) of bearings 2 and 3 are rotatinq in 

stagnant fluid. 


3) Fluid is rotating at cage speed between bearings 1 and 2 and between bear- 
ings 3 and 4. 

4) The power dissipated by the shaft section between bearings 2 and 3 was 

determined by averaging the power estimates based on inlet fluid properties 
and exit fluid properties. ■ . 

5) The power estimated for disks and cylinders can be superimposed to obtain 
the combined slmger-shaft power. 


bpaH ?* "I 1 S 1Qure 4 ’ 3 ' 6 ’ onl -y the bearing frictional power for the loaded 
Fn!^w 9 v i * 3 ' + W u 3S vaned as a faction of bearing temperature differences, 
hnergy values for other components downstream of bearings 2 and 3 would vary as 
additional heat is added to the coolant. This effect was considered secondary in 
^° ther unc ® rta ?"ties in the analysis. The other side of the energy 
d “ ls the e ^rgy absorbed by the bearing coolant flow. This was estimated, and 
the summary results are provided in Figure 4.3.5. As shown, 182.1 kW of energy was 

t0 5 ab ^ orbed b y the bearing coolant flow. Since this is in excess of the 
esbima b ed tG be generated, it is suggested that a temperature gradient 
across the bearings could explain the discrepancy. Figure 4.3.7 illustrates, in 
graphical form, the information provided in Figure 4.3.6. Also shown in Figure 4 3 7 

is the temperature gradient across the bearings (2 and 3) necessary to achieve an 
overall energy balance. 
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In summary, the results of the analysis strongly indicate that the inner races 
of bearings 2 and 3 are considerably higher in temperature than the outer races. 
This condition can become unstable causing excessive internal loads and temperatures 
resulting in rapid deterioration and failure of the bearing. In addition, careful 
consideration must be taken to ensure the unloaded bearings (1 and 4) do not become 
unloaded to the extent that significant slip or skidding can occur between the balls 
and inner races. 

4.3.3 Estimation of Heat Loads for the BMT Thermal Model 


Heat is generated in the BMT from mechanical friction and viscous work done on 
the coolant. To realistically estimate bearing component temperature, this energy 
must be evaluated and applied correctly in the BMT thermal model. An estimate of the 
total energy generated can be obtained from an overall energy balance on previous 
test data. This was done in Section 4.3.2and is shown in Figure 4.3.5. Once the 
total energy is determined, the problem remaining is to correctly evaluate the 
distribution of viscous and frictional energy. The bearing contact surface tempera- 
tures are obviously much more sensitive to frictional energy changes than changes in 
viscous energy. 


The distribution of energy for the BMT thermal model, with LFL coolant, was 
evaluated as follows: L 

1) Estimate total energy generated from a heat balance on the coolant (see 
Figure 4.3.5 of Section 4.3.2). 

2) Calculate the viscous energy generated by each component in the tester (see 
Figure 4.3.6 of Section 4.3.2). 

3) Using measured torque values, estimate appropriate friction coefficients, 
and calculate friction energy distribution (i.e. inner race, outer race). 

4) Adjust viscous energy distribution such that the total input matches the 
overall heat balance determined from test data. 

Steps 1 and 2 had already been done and the results provided in Section 4.3.2. 

Battel le (Reference 34) measured the bearing torque as a function of axial load 
and surface conditions. Figure 4.3.8 is a reproduction of their data. Also shown, 
for comparison, is the bearing torque predicted by SHABERTH I. The predicted torque 
curve was generated using a friction coefficient of 0.2. The predicted curve is rea- 
sonably close to the PTFE burnished condition. The referenced value of friction 
coefficient for this surface condition is .08. The reason for this difference is not 
evident, however, since power and heat rate are equivalent, a coefficient of friction 
of 0.2 in the SHABERTH program provides a value reasonably close to the measurement. 
A further verification of the analysis is shown in Figure 4.3.9. This is data taken 
from the same reference, compared with results from the SHABERTH I bearing analysis 
computer code. Although the calculated and measured axial deflections diverge with 
increased load, the contact angles agree fairly well. Consequently a friction 
coefficient of 0.2 was used in SHABERTH to obtain frictional heat generation values 
for the thermal model. 

To complete the frictional heat evaluation the bearing load as well as the fric- 
tion coefficient must be known. There is uncertainty in the axial load seen by the 
bearings due to the unknown pressure loss across the carriers. According to the 
measured AP the axial load on the bearing should be about 2500 pounds. Although this 
value is suspect because of fluid vortices, it is probably the most representative at 
this time. 
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FIGURE 4.3.8 BEARING TORQUE AS A FUNCTION OF AXIAL LOAD 




FIGURE 4.3.9 BEARING DEFLECTIONS AND CONTACT ANGLE AS A FUNCTION OF AXIAL LOAD 
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Therefore this axial load in combination with a friction value of 0.2 was used 
to estimate the frictional heat generated. The amount of friction energy was sub- 
tracted from the total measured energy, and the remainder was assumed to be viscous. 
The viscous energy was distributed in the proper proportions as indicated by the 
viscous work analysis done in Section 4.3.2. It was assumed that no lubrication was 
provided by the LN- coolant. This assumption will be reevaluated for the analyses of 
the LOX coolant. 

Another important parameter, with a high degree of uncertainty, is the two phase 
heat transfer coefficient. There is no experimental data on this parameter that 
corresponds to the operating conditions of the bearing. Therefore, extrapolations 
have been made to estimate these parameters for the thermal model. The uncertainty in 
these extrapolations are unknown. 

The above points out two critical areas of cryogenic bearing analysis requiring 
experimental work: 1) Evaluation of the friction coefficient for high speed bear- 
ings using film transfer lubrication in cryogenics, and 2) Evaluation of two phase 
heat transfer coefficients for high speed bearings in cryogenics. 

4.3.4 Iterative Determination of Distributed Frictional Heat Generation in Bearing 
Contacts 


The bearing thermal model has the capability of allowing the frictional heat 
generated in the rolling element/race contact to be distributed along the major axis 
of the contact ellipse. The SHABERTH bearing analysis code provides the frictional 
heat generated in the bearing contacts. The calculation technique used by the code 
solves for the friction force and slip velocity at finite increments in the contact 
ellipse. The sum of the product of these parameters is the friction heat generated 
in the contact. By modification of the SHABERTH code, the heat generation rate 
distribution along the major axis of the contact ellipse can be obtained. Examples 
of these distributions are shown in Figures 4.3.10 and 4.3.11 for the inner and outer 
races respectively. It is of interest to observe that the maximum normal force does 
not coincide with the location of maximum heat generation. The most heavily stressed 
area of the contact is therefore not at the highest contact temperature. This is a 
fortunate condition for those cases where material properties become marginal due to 
high temperatures. 

An iteration procedure has been developed to determine the temperature gradient 
across the BMT bearing utilizing the "SHABERTH" bearing program and the Bearing 
Thermal Model with "SINDA" used for analyzing the thermal system represented in 
electrical analog or lumped parameters form. A specific axial reaction is selected 
for analysis and put into the "SHABERTH" data base which is set up with configuration 
data for a 57 mm bearing. A first estimate temperature gradient across the bearing 
is assumed and added to the data base along with appropriate temperature varying 
thermal property data for the shaft, inner race, rolling element, outer race and 
housing. A friction coefficient of 0.2 was used. The "SHABERTH" model is then 
exercised with this data base and the resultant information contains the frictional 
heat generated at the track areas of the races and rolling element for a specific 
axial reaction. This heat, generated from mechanical friction, is input into the 
Bearing Thermal Model data base along with heat generated from viscous work done on 
the coolant. 

The thermal model is then exercised for a specific data base with the resulting 
output being steady state temperatures for each node. The average temperature for 
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each bearing component (i.e., shaft, inner race, rolling element, outer race, and 
housing is then calculated as shown 



i=N 



T. = Temp, of Node 
vj = Volume of Node 
N 1 = Total Number of Nodes 
Vy = Total Volume of N Nodes 


The Bearing Thermal Model has been coded to calculate the average temperatures 
of the bearing components and to determine the thermal property data needed from 
derived curve fit equations. 


The average temperature values for the bearing components are then placed back 
into the "SHABERTH" data base and replace the original assumed values. Appropriate 
temperature varying thermal property data is also replaced in the data base. The 
"SHABERTH" program is then executed again to determine new distributed frictional 
heat generation data. This data is compared with the previously obtained values to 
determine whether another iteration is required. It is important to obtain not only 
the correct AT across the bearing, but also the correct bearing component tempera- 
tures. The aT across the bearing changes the bearing internal clearance. Reduced 
clearance results in a larger load on the bearing which in turn increases the fric- 
tional heat generated. The individual component temperatures are used to determine 
such thermal property data as the coefficient of thermal expansion and the modulus of 
elasticity, the values of which, also affect the bearing internal clearances. 

The iteration procedure outlined above was performed for an axial reaction of 
2500 pounds. The operating parameters; inlet coolant temperature of -176°C and mass 
flow rate of 3.2 LB/SEC per bearing set, were selected for a "worst case" situation. 
The results of this iteration are shown in Figure 4.3.12. The preliminary frictional 
heat generation data was obtained from a "SHABERTH" run with an initial assumed aT 
across the bearing of 200 C. This temperature gradient was then resubmitted into the 
"SHABERTH" bearing data base along with a friction coefficient of 0.2. The first 
iteration frictional heat generation data is significantly smaller than the prelimi- 
nary data. Therefore, the Bearing Thermal Model does not support the original 
assumption of a 200 C aT across the bearing. A second iteration shows the heat 
generation data holding relatively constant which supports a preliminary AT across 
the bearing of 25 C and average temperatures of -137°C, -122°C and -162°C for the 
inner race, rolling element, and outer race respectively. This indicates a stable 
set of thermal conditions, i.e., the heat generated and bearing temperatures do not 
increase without limit. 


4.4 Heat Transfer Coefficients for Bearing Thermal Model 
4.4.1 Two-Phase Heat Transfer in Cryogenic Systems 

In any attempt to predict failure using the bearing thermal model, the effects 
of surface boiling on the coolant performance must be considered. Preliminary SINDA 
runs showed, as expected, several 'hot spots' near the ball -race contact areas where 
frictional heat generation is high. Local boiling can occur here (even though the 
bulk of the fluid is subcooled) causing an abrupt reduction in the surface heat 
transfer coefficient and potential increases over the maximum heat flux temperature 
difference. Therefore, a boiling curve (heat transfer coefficient versus surface 
temperature) must be established. 
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Test data or published correlations for heat transfer coefficients for sub- 
cooled, high pressure, high flow, boiling cryogenic systems have not been found in 
the literature (See Reference 16-18, 20-33, and 38-40). There appears to be no test 
data that matches the thermal and flow conditions in the bearing tester. Therefore, 
available information must be extrapolated to the bearing tester conditions and later 
verified by data from the tester. 

The following describes the rationale for extrapolating existing data/correla- 
tion to represent the thermal and flow conditions in the bearing tester. Although 
similar techniques will be used for liquid oxygen, the results of the following 
development are for LN,,. 

The forced convection subcooled boiling process can be separated into four 
phases: 

1) Forced convection to single liquid phase. 

2) A combination of forced convection and nucleate boiling until a maximum 
heat flux is reached. This maximum heat flux is commonly referred to as 
the "burnout" heat flux. 

3} Partial film boiling. 

4) Stable film boiling. 

The forced convection phase occurs without vapor generation and this phase of the 
heat transfer process is fairly well represented with published correlations and test 
data. 


Nucleate boiling is characterized by a very high heat transfer rate for small 
temperature differences. There are two subregimes in nucleate boiling: local 
boiling and bulk boiling. Local boiling is nucleate boiling in a subcooled liquid 
which is characteristic of the conditions in the bearing tester. Bulk boiling is 
nucleate boiling in a saturated liquid which will not occur in the tester unless the 
outlet fluid temperature reaches the saturated value. As will be discussed later, 
the effects of local nucleate boiling are significantly influenced by flow velocity 
and subcooling. When the population of bubbles becomes too high the bubbles leaving 
the surface restrict the flow of liquid to the surface and insulate it thus causing a 
rapid rise in surface temperature. This is defined as the "burnout" heat flux. It 
is important to define this point because in general the heat transfer mechanism 
changes from nucleate boiling to stable film boiling with a significant increase in 
surface temperature. Stable operation in the partial film boiling regime generally 
does not occur. Generally when the peak flux or "burnout" point has been passed the 
system moves to the stable film boiling regime. 

In order to account for these various heat transfer mechanisms in the bearing 
thermal model, criteria was established for the changing modes of heat transfer as 
functions of the surface temperature of the specific nodal network. For example, as 
the inner race or ball surface node temperature increase from the subcooled liquid 
temperature to temperatures considerably in excess of the fluid saturation 
temperature, all the above modes of heat transfer will have been encountered. 
Consequently, the node surface temperature range for each mode of heat transfer must 
be established and the appropriate heat transfer correlation determined. 

FORCED CONVECTION HEAT TRANSFER 

As long as the surface temperature is below the local fluid saturation tempera- 
ture, forced convection is the heat transfer mechanism. The appropriate heat trans- 
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fer correlation for these conditions is the Dittus - Boelter equation for the inner 
and outer race surfaces and other surfaces in the flow stream with the exception of 
the ball. The heat transfer correlation for the ball was taken from Reference 40 
which correlates a wide range of heat transfer data for spheres. 


The Dittus - Boelter equation is as follows:* 

D h c D e = .023 N re ' 8 N pr ' 4 
k 

And Katsnel Ison's equation for spheres (Ref 40) is: 


2) ^ = 2 + 0.03 N or 0 - 33 N re °‘ 54 + 0.35N °- 356 N re 0<58 

k v 


FORCED CONVECTION NUCLEATE BOILING 


In this heat transfer regime, both forced convection and subcooled boiling occur 
simultaneously and no generalized correlations, for subcooled cryogenics, have been 
found. The first steo in developing an approach for extropolation of available 
information to these conditions is to select a pool boiling correlation that matches 
test data for high pressure LN ? . There are several pool boiling correlations avail- 
able to select from. Examples are: Rohsenow, Forster - Zuber, Borishanskiy - 
Minchenko, Kutateladze and others. Reference 38 compares these correlations with 
nucleate boiling test data for several cryogens including LN^. Figure 4.4.1 shows a 
sample of those comparisons. As shown, there is a considerable pressure effect on 
the boiling heat flux at a given value of AT. Also shown is the superior correlation 
with test data of the Russian equation compared to Rohsenhow's expression. Based on 
these results, the Borishanskiy - Minchenko correlation was selected as the best 
representation of pool boiling LN^ data. This equation can be expressed as follows: 

3) T„ - T s - f(P) q b °' 3 


The next step in defining the heat transfer regimes is to determine the node 
surface temperature at the peak heat flux or "burnout" point. An expression, devel- 
oped by Kutateladze, was selected based on the comparison with pool boiling burnout 
flux for LN^ provided in Reference 38. This expression is as follows: 

4) °- 16 f(p) 

Since the bearing tester operates in subcooled liquid, the effects of subcooling 
on heat flux should be accounted for. Reference 39 provides the following expression 
for correcting saturated pool boiling heat flux for subcooling. 

5) %b = 1 + 0 [T s - T L ] 

% 

The above boiling data correlations are for static boiling systems, wherein the 
flow velocity is small or non existent. The procedure for extrapolating these corre- 
lations to a high flow subcooled boiling cryogenic system is as follows. 

* Refer to nomenclature on page 59 for definition of terms. 
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1) Estimate the surface temperature at which nucleate boiling initiates in the flow 
system. 

2) Estimate the surface temperature at the burnout or maximum heat flux conditions. 
After the above conditions are established it will be assumed that: 

1) Forced convection is the dominant mode of heat transfer prior to the surface 
temperature reaching the incipient nucleate boiling temperature. 

2) When the surface temperature is between the incipient boiling and the tempera- 
ture corresponding to peak boiling, a combination of forced convection and 
nucleate boiling occurs. 

3) When the surface temperature exceeds the temperature corresponding to the peak 
flux condition, stable film boiling occurs and this condition can be represented 
by the forced convection relationships using vapor properties. 

ESTIMATION OF SURFACE TEMPERATURE FOR BOILING INITIATION 

It is assumed that nucleate boiling will begin when the boiling heat flux as 
represented by equation 3 equals the heat that can be removed by forced convection. 
For example: 

6 > - h c<‘w - T L> 


Then from equation 3: 

7) T w - T s * f ( p > thcOw - t l>1 °' 3 

In the above equation Tw is the surface temperature required to initiate 
nucleate boiling. 

ESTIMATION OF SURFACE TEMPERATURE AT MAXIMUM HEATING CONDITIONS 

The pool boiling "burnout" heat flux is estimated from equation 4 and corrected 
for the appropriate subcooling by equation 5. For the high flow system both boiling 
and forced convection occur throughout the nucleate boiling regime. Therefore, at 
burnout the total heat flux is assumed to be a combination of these mechanisms. The 
following expression can be written assuming that the nucleate boiling relationship 
still holds. 

8) W s .f<p) i h c< T w * V + <W 0 ' 3 


Tw in this case is the wall temperature at "burnout" and when this temperature 
is exceeded stable film boiling is assumed to occur. 

The above technique has been programmed for computer evaluation of the thermal 
conductors required in the SINDA bearing thermal model. In general, it was found 
that forced convection heat transfer dominated in the high velocity regions. Heat 
transfer from the high-speed components, inner races, balls, etc, is affected very 
little by nucleate boiling. The major effect is at the point where the surface is 
estimated to be vapor blanketed (i.e., the point of "burnout"). 
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NOMENCLATURE FOR SECTION 4.4.1 
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The magnitude of the "burnout" heat flux can also be significantly affected by 
surface conditions such as roughness, coatings, etc. The maximum heat flux condition 
can be reached prematurely if flow instabilities or oscillation occurs. These rapid 
changes in flow can lower the local pressure increasing the local superheat near the 
wall and cause rapidly increased boiling. These added uncertainties along with the 
lack of experimental data relative to the thermodynamic conditions in the bearing 
tester emphasizes the importance of experimental data from the tester. By comparing 
measured and calculated temperature data the validity of the thermal model can be as- 
sessed. 


SENSITIVITY OF BEARING COMPONENT TEMPERATURES TO TRANSITION TO FILM BOILING 

Temperature measurements can serve as a guide in evaluating the transition point 
from forced convection nucleate- boiling (burnout flux) to film boiling. An example 
of the importance in evaluating this parameter is shown in Figure 4.4.2. The average 
bearing component temperatures are shown as a function of the maximum wall superheat 
(Tw-Ts) MA w for forced convection nucleate boiling. As the maximum superheat is in- 
creased, Considerably more surface area of the bearing components remains in the 
nucleate boiling regime. Therefore, these areas do not become vapor blanketed and 
the surface-to-fluid heat transfer remains relatively high resulting in lower temper- 
atures. As. ( Tw-T s )mav is increased to the point that all surfaces remain in the 
nucleate boiling regime, the component temperatures become independent of the transi- 
tion from nucleate to film boiling. In other words, for a given heat load, the 
transition point can be increased such that the surface temperature does not reach 
the transition point and film boiling does not occur. In Figure 4.4.2 this is 
represented by the higher AT MA) , values where the bearing component temperatures 
become independent of aT^. a 

The analysis described herein predicts that the peak flux will be reached at 
surface temperatures between 2.6 and 3.3°R, above saturation, depending on the 
magnitude of the forced convection heat transfer coefficient. Also shown is a 
temperature measurement for the number 2 bearing outer race. This measurement is in 
reasonable agreement with the calculated outer race temperature at the estimated 
4 Tmax point. 

4.4.2 Fluid Heat Transfer Regimes and Film Coefficie nts for LOX and LN„ 

Techniques developed in Section 4.4,1 were used to evaluate the surface tempera- 
ture for the incipient nucleate boiling range for LOX. The heat transfer regimes as 
a function of wall superheat (T -T ) are shown in Figure 4.4.3. As shown, nucleate 
boiling for the forced convection COX flow begins for a wall superheat of 12.8°R and 
the. "burnout" heat flux occurs at a wall superheat of 13.8°R. This indicates that 
typical bearing urfaces will become vapor blanketed when the surface temperature 
exceeds the LOX saturation temperature by about 13°R, This compares with a value of 
about 3°R for LN£ coolant as shown in Figure 4.4.4. 

4.4.3 Transition to Film Boiling Temperature for LN,, 

An analysis was performed using the Bearing Thermal Model to evaluate the 
transition temperature at which LN~ film boiling starts in the BMT bearing component. 
When the transition superheat (3°F7 predicted by an earlier calculation was used, the 
steady state outer race temperatures compared favorably with the limited BMT thermo- 
couple data. Therefore, this value (&T = 3°F) will continue to be used until more 
test data becomes available. The bearing analysis assumed an axial load of 2,500 
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FIGURE <.4.1 COMPARISON OF NUCLEATE BOILING CORRELATIONS 
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lbs, an LN ? flow rate of 6.3 lbs/sec brg set, with an inlet temperature of -305°F. 
The stead/ state average component temperatures (bearing #2) at the transition 
superheats were: 


Ball -185°F 
Outer Race -252°F 
Inner Race -262°F 
Shaft -272°F. 


Although the favorable agreement of analysis and test data is encouraging, 
considerable additional test data at varied operating conditions are required for 
bearing model verifications. 


The Bearing Thermal Model was exercised for various transition temperatures 
using LOX as the coolant. The transition temperature for the bearing - coolant 
combination is that temperature at which the maximum surface boiling heat flux 
occurs. This temperature is also called the "burnout" temperature because a higher 
heat flux causes the surface to dry out causing it to become blanketed with a vapor 
film and the surface temperature to drastically increase. This transition tempera- 
ture can be expressed as transition superheat - the difference between the surface 
temperature and the fluid's saturation temperature at the maximum boiling heat flux.. 

Analyses results are shown in Figures 4.4.5 through 4.4.7. Figure 4.4.5 shows 
the average, steady state component temperatures as a function of the transition 
superheat. The average temperatures as a function of the transition superheat. The 
average temperatures are not very sensitive to transition superheat changes in the 
range investigated (aT = 2, 15°F). Figure 4.4.4 shows the temperatures of the two 
outer race nodes whose location most closely corresponds to the placement of the 
thermocouples in the BMT. Also shown, for comparison, is an outer race temperature 
measured on the BMT under similar conditions. A previous calculation (Section 
4.4.1), predicted a transition superheat of about 3°F for LN„. This value gives a 
very close match between node #463 and the measured temperature of -268°F. There- 
fore, we will continue to use this value (AT - 3°F) for the transition superheat in 
future studies until more test data is available and a more comprehensive correlation 
can be established. Figure 4.4.5 shows the maximum track node temperatures for the 
races and the ball. At the transition superheat of 3°, all the track nodes have 
surface temperatures well above the LN,, saturation temperature of -241°F. Conse- 
quently, the loaded contact surfaces are' 1 vapor blanketed and receive little, if any, 
lubrication from the LN^. At the transition superheat of 3°F, the track nodes had 
the following temperatures: 


Ball Track +130°F 
Inner Race Track +52°F 
Outer Rack Track -132°F 


The bearing pair analyzed has an applied axial load of 2,500 lbs, an flow 
rate of 6.3 lb/sec brg set, and an inlet temperature of -305°F. 

4.4.4 Transition to Film Boiling Temperature for LOX 

The transition temperature for the bearing - LOX combination is that temperature 
at which the maximum forced convection nucleate boiling heat flux occurs and subse- 
quently causes stable film boiling to take place on the bearing surfaces. This can 
also be expressed as a temperature difference between the surface and the saturation 
temperature of the fluid (superheat) at the maximum heat flux. 


63 



spectra research systems 


ORIGINAL PAGE 
OF POOR QUALtTY 




IICIIRI I A, A COIPOHEHT TEIIPEMIUBE AS A FUIOIOK OF I'ACL SUPERHEAT 
AT IIAI FEUI CONOITtOMS 


GBF. 12 



5 1(1 15 20 T5 


(I IJAU ' 'sAT 1 AT m * "W riu >- * r 


r iraj " E M.9 CBlPOXroi TBIPERAIUflE AS A FURCTtOA or HALL SUPERHEAT 
AT IIAl FEW COnSITIOHS 


FI GUBC A.A.10 r RACK NODE TEMPERATURES AS A FlaCTlOK OF UAEL SUPERHEAT 
AT RAT FUJI C0TI0IT10NS 



64 


spectra research systems 


. ^ second parameter', the temperature at which the vapor properties were evaluated 
in the fluid heat transfer coefficient calculations, was investigated. Durinq vapor 
generation on a surface, the temperature gradient between the surface and the 
subcooled liquid is unknown, and the vapor properties must be evaluated at some 
reference temperature. Two possibilities were examined: (1) A film temperature 
equal to the average of the surface temperature and the fluid's saturation tempera- 
ture, and (2) The saturation temperature. 

The analysis was performed using the Bearing Thermal Model to evaluate the 
temperature at which the transition to film boiling occurs for LQX in the BMT. Vapor 
properties needed for the conductance calculations were evaluated at both film 
temperature and saturation temperature. The resulting outer race temperatures were 
compared to the BMT thermocouple data for the LOX test to provide a comparison of 
predicted and measured temperatures. A superheat of 5°F with the vapor properties 
evaluated at the film temperature was chosen for the transition temperature. This 
produced an outer race temperature for Bearing § 2 of -217°F and the following average 
component temperatures: y y 


Ball 

Outer Race 
Inner Race 
Shaft 


- 96°F 
-207°F 
-200°F 
-227°F. 


This bearing pair had an applied axial load of 2500 lb, a LOX flow rate of 4.5 Ib/sec 
and a LOX inlet temperature of -268°F. 7 


The results of this effort provide detailed thermal and mechanical 
support the design and performance evaluation of high-speed angular contact ball 
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properties produces a strong difference in these curves. Using the chosen superheat 
of 5 C F at burnout, the track nodes have the following steady state temperatures : 

Ball Track +257° F 

Inner Race Track +170° F 

Outer Race Track -10° F. 

While these temperatures are the highest in the contact areas, all the track node 
temperatures are well above the saturation temperature of the LOX. Consequently the 
loaded contact surfaces are vapor blanketed and receive very little if any film 
lubrication from the LOX. This bearing pair had a 2500 lb applied axial load, a LOX 
flow rate of 4.5 lb/sec, and an inlet coolant temperature of -268° F. 

4. 5 Software Modifications and Development 


Two large computer programs, SHABERTH and SINDA were used for many of the 
analyses conducted during the contract period. SINDA, a numerical analysis program 
designed specifically for solving resistor-capacitor network representations of 
thermal systems by forward -backward differencing methods, was used to build the SRS 
Bearing Thermal Model data base. Appendix A contains a representative data base and 
illustrates the variety and magnitude of data considered in the thermal model. 
Appendix B illustrates the basic steady-state output temperature data. In addition 
to the basic node data and conductor data common to most thermal analyses, the SRS 
Bearing Thermal Model data base contains temperature varying thermal property data, 
capacitance values for each node to allow transient analyses, time varying heat 
source capabilities, as well as data manipulation capabilities. SHABERTH, a large 
program that predicts general bearing performance characteristics, was used to model 
the Bearing and Materials Tester using geometrical dimensions from the original LOX 
tester. A representative input data file and output data file for this program is 
shown in Appendix C. The SHABERTH program considers such factors of bearing opera- 
tion as clearances, pitch diameter, rail-land width and diameter, cage weight, number 
of rolling elements, number of bearings, steel data, ball diameter, curvatures, CLA 
roughness, RMS. asperity slope, friction coefficients, fit data, axial and radial 
loading, material properties, temperature gradients across bearing, shaft geometry 
and is considered a state-of-the-art bearing analysis program. SHABERTH output 
contains such information as axial and radial shaft deflections, reaction forces, 
moments, frictional heat generation, friction torque, bearing operating clearances, 
cage speed data, rolling element angular speeds, speed vector angles, normal fo.rces. 
Hertz stresses, and contact angles. 

Many modifications were made to these programs to adapt them to the specific 
problems of the BMT. Also, several other smaller routines were developed to enhance 
the capabilities of both SHABERTH and SINDA. When SHABERTH was first acquired, it 
would not run under the current version of Fortran available on the MSFC Uni vac 
computer. The code was finally updated to Fortran V and is now fully updated to 
ASCII Fortran (also known as Fortran 77). Besides the code changes necessary to 
update the SHABERTH code, two subroutines, FRINT and FMIX, were modified to store the 
characterise cs of the ball -race contact ellipse. SHABERTH models these ellipses by 
slices. The width of each slice, the frictional heat generated at each slice and the 
major and minor axis dimension of both the inner and outer race contact ellipse are 
stored for use in the SRS Bearing Thermal Model. 

SINDA is capable of computing both steady-state and transient heat source 
temperature profiles and allows for the addition of user-written subroutines to be 
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added and compiled along with the standard library of subroutines available. A 
routine has been added to compute the mechanical and thermal properties of the 
bearing set at the predicted temperatures. The subroutine, called VARIABLES 2, reads 
a file that contains the volumes of each node in the SRS Bearing Thermal Model. It 
then computes the average temperature of each component of the bearing (inner race, 
outer race, and rolling element) based on these volumes. The modulus of elasticity 
and coefficient of thermal expansion for each bearing component is then calculated. 
VARIABLES 2 is included in the SINDA input file listed in Appendix A. 

The SRS Bearing Thermal Model was designed to be adaptable to changing various 
operating parameters such as coolant flow rate, coolant inlet temperature, and axial 
load for either LOX or LN„ coolant. However, these changes require re-calculation of 
over 1,000 conductance values for each run. To calculate each by hand would be a 
laborious and time consuming process. The following special routines were developed 
to facilitate this process and are listed in Appendix D: COND/SRS, HEAT/SRS* SRS- 
DATA, and UPDATE/SRS. These subroutines have been combined with SHABERTH to form a 
new program called SHABERTH/SRS that updates a Bearing Thermal Model file quick \y and 
efficiently when new operating conditions are being considered. The rootlet and 
their functions are described below. 

COND/SRS - This subroutine calculates the surface-to-fluid conductances for the SRS 
Bearing Thermal Model. The heat transfer coefficient is also calculated using 
methods outlined in Section 4.4. For each node, a table of conductances vs, ‘smpera- 
ture is produced that accounts for the changing conductances through the different 
heat transfer phases (forced convection, boiling, and film boiling). The output of 
COND/SRS is a list of these tables arranged in SINDA input format. SINDA can then 
interpolate through these tables to get the conductance for the surface nodes as a 
function of temperature. C0ND/PR0G also adds commends to the SINDA file that identi- 
fy the run and describe the operating conditions. 

HEAT/ SRS - This subroutine distributes the frictional heat generated at the bearing 
contact surfaces, that is calculated by SHABERTH, to the proper nodes on the ball and 
races in the SRS Thermal Model data base. As mentioned above, the FRINT and FMIX 
subroutines of SHABERTH isolate the frictional heat generated into slices of the 
contact ellipse. HEAT/SRS also calculates the metal -to-metal conductances between 
the ball and the races since the contact surface area is dependent upon the load 
applied. 


UPDATE/SRS - This subroutine takes the output of COND/SRS and HEAT/SRS and updates a 
SRS Bearing Thermal Model input file with the re-calculated conductances and fric- 
tional heat generation. 


SRSDATA - This is a FORTRAN BLGCK DATA subroutine that contains all the constants 
needed by the COND/SRS and HEAT/SRS, including various LOX and LN„ properties (densi- 
ty, specific heat, thermal conductivity, Prandtl number, viscosity), the thermal 
conductivities for the four metals used in the BMT (Inconel, Armalon, 440-c Stainless 
Steel, Waspaloy), and several node properties of the Bearing Thermal Model (surface 
area, depth, location relative to shaft center, hydraulic diameter of the passage at 
the node's surface, and a relative velocity indicator). This data is held in COMMON 
blocks that the other subroutines can access. 
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5.0 BEARING ANALYSES 


This section provides detailed descriptions of specific tasks accomplished and 
the results obtained from the bearing modeling efforts. These analyses were 
performed in support of the Bearing and Materials Tester Program, as required, and do 
not represent the total analysis effort for the design and development of the BMT. 

5. 1 Bearing Life Prediction Analysis 

5.1.1 Modific ation of Bearing Life Formulas for Rolling Bearings Operating in 
Cryogeni c Envi ronment ~ ~ ~ ' “ 

Accepted techniques (4,6,12) have been developed to predict life for well 
lubricated bearings. These techniques lack experimental verification for high-speed 
rolling bearings operating in cryogenics, it is therefore advantageous to investigate 
the potential for modifying current bearing life formulas for rolling bearings 
operating in a cryogenic environment. Typical cryogenics (LH ? , LN ? LOX) have good 
heat removal characteristics; however, they lack the viscous properries necessary to 
provide a fluid film of sufficient strength to prevent metal -to-metal contact between 
the rolling elements and the race. The degree these films support the interface 
loads has not been experimentally determined. 

Another uncertainty is the role cryogen film lubrication plays when solid lubri- 
cants are used to provide a lubricating surface film. The interaction between the 
solid film and fluid film are not known. This can be further complicated when the 
contact surface temperature exceeds the saturation temperature of the coolant to the 
extent that the surfaces become vapor blanketed and the surface film is that of a 
gas. Some experiments (15) indicate that LOX is potentially a better lubricant than 
non-reactive cryogens such as LHp and LNp. This results from the formation of 
surface films (oxides) that provide a measures of protection for the metal surfaces in 
contact. Here again, the rate of oxide formation can be a function of surface 
temperature and the degree of protection can be a function of the sliding velocity or 
rate of film removal. Considerably more work is needed to establish the film lubrica- 
tion properties of various bearing materials and solid lubricants operating in 
cryogenics. 


The classical failure mode of a rolling contact between hard steel surfaces is 
contact, fatigue. The life formula as developed by Lundberg (4) is based -on the 
assumption that the failure is initiated by a decisive stress amplitude located at 
some distance below the surface. The relationship between loads, stresses, and 
stress location is developed from the bearing dynamics and Hertz theory of contact 
stresses. Bearing life dispersion has been obtained by applying the Weibull dis- 
tribution function to bearing life data. Experience has shown that this technique 
adequately predicts bearing life for well lubricated bearings. Since better lubri- 
cants and materials have become available, the original life formula has been mod- 
ified by factors to account for the improved life derived from these improvements 
(6,12). These factors, however, do not account for the conditions encountered in 
cryogenic operation. As an example, the development of stresses from the Hertz 
theory does not account for tangential forces. Since the contact friction forces 
with cryogenic fluids is expected to be significantly larger than the forces encoun- 
tered in oil lubricated bearings, the maximum shear force will probably be at, or 
near, the surface rather than below the surface. Furthermore, the solid film lubri- 
cant will enhance the asperity friction. Consequently the apparent contact friction 
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force should fall somewhere between dry asperity friction and hydrodynamic film 
friction. This presents a problem in determining the life adjustment factor for 
lubrication that will account for the degradation of lubrication capability of the 
cryogenic. 

The SHABERTH I bearing analysis computer program has been used as an aid in 
predicting bearing life. This program uses the life formula developed in reference 4 
with life adjustment factors to account for improved materials and lubrication film 
thickness. The life adjustment factors will be adjusted to represent cryogenic 
operations as data becomes available from the bearing tester program. 

SHABERTH I uses the following technique for estimating bearing life: 



^IQm 1S racewa y fatigue life in millions of revolutions, as determined by 
Lundberg and Palmgren (4), and Q is the dynamic capacity which is the load for 
which the bearing raceway will have 90% assurance of surviving one million revo- 
lutions. Q is the raceway equivalent load. Q and Q are calculated by the 
SHABERTH colter program. cm em 

Since experience has shown, for ball bearings, that the races usually fail 
before the balls due to the changing axis of rotation as the balls orbit the shaft, 
the bearing fatigue life is based on the life of the races. The L. n life for a 
raceway is given by: iU 

L 10n = a 2 a 3 a 3* L 10m where 

a 2 is a life improvement factor to account for improved materials, 

a 3 is a life improvement factor to account for full film lubrication, and 

a 3* 1S a f ac t° r ^ ess than one when full film lubrication is not obtained. 

a 3 * is a reduction factor accounting for the effect of surface asperity inter- 
action and was deduced from a best fit to a large body of rolling contact life test 
data (12). This function was extracted from the SHABERTH I program and is shown in 
Figure 5.1.1. Also shown is the lubricant life factor used by SHABERTH II, which 
will be discussed later. As shown, the SHABERTH I lubrication factor is always <1 
and is therefore never a life improvement factor. If the film thickness to surface 
roughness ration (h/a) is greater than about 3.5, a life improvement factor >1 can be 
input (a 2 ). 

A limiting value of about 0.52 is shown for h/o approaching zero, or zero film 
thickness. Since this condition is believed to be approached in LN„, LOX, and LH 0 , 
this is the value that would be applied with these fluids. The detemi nation of this 
life factor did not, however, consider the benefits of the solid lubricant film 
transfer, system used to lubricate the tester bearings. The effectiveness of this 
lubrication on bearing life can be evaluated from the bearing and materials tester 
data. SHABERTH II has, in effect, combined the life reduction and improvement due to 
lubrication into one function. This function is developed in Reference 6. The 
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FIGURE 5.1.1 LUBRICATION LIFE ADJUSTMENT FACTORS 
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important point is that for b/a approaching zero, the life reduction factor is 0.21, 
which is less than half the factor used in SHABERTH I. This, again, illustrates the 
need for experimental data for cryogenic operation in these low film regions. 

is the raceway fatigue life for 90% survivability, 

ie I_iq life for thp hparirm ic 
1/e 


is: 


As stated earlier, the l 
a the 

Eq 2 


Since the bearing has two race¥Hhe I_ 10 life for the bearing is 


•10b 


-e 

L 10l * L 102 


"1 

02 J 


where 1, 2 denote inner and outer races. 

Since one bearing set consists of two bearings, the L 1Q life of the bearing set 

Eq 3 r 1 " 1/S 


-10bs 
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L 10bl + L 10b2 
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And the life of the total bearing system is; 
Eq 4 r _ _ n "1/e 


'lObst 




-e 


-e 


L 10bsi + L 10bs2 


. Equati on 4 is the 90% survivabil ity life for all the bearings on the tester 
shaft. Since a 10* failure rate is not acceptable for Shuttle components. Equation 4 
^Li? 4 - m0d * lf i ed t0 predict probabilities of success greater than 99%. Techniques to 
correlate failure rates in these low regimes are addressed in References 1 and 5. 
These techniques suffer from the same deficiencies as discussed earlier when applied 
to bearings operating in low viscosity fluids such as LN os LOX and LH 0 . However, 

they do provide a technique for evaluating the data derived from the bearing and 
materials Tester. 

The previous discussion has addressed fatigue type failures which is the domi- 
nant failure mode f or we11 designed and lubricated bearings. Obviously, there are 
other type failures such as excessive wear, overheat cracking, overload crackinq, de- 
struction by material softening due to unstable overheating, etc. These type fail- 
ures must be overcome before life data relating to bearing fatigue can be obtained. 


to 


The effort to date has provided a bibliography of reference material pertaininq 
bearing fatigue life (References 3 through 15). Life data relating to bearinqs 
operating in cryogenics is scarce and relating existinq semi -empirical life formulas 

dl i 1 f s . re< i u J re additional test data. The Weibull distribution does not 
represent the data in the region of reliabilities >90%. Techniques have been found 

™ r * ent 1n these re 9 imes and should, with modification, represent 
the data obtained from the cryogenic tests. 

5 «1- 2 Cumulative Load Assessment for BMT Bearing Life 

A basic problem in fatigue analysis is the prediction of life for an element 
when time-dependent loading is applied. By progressively loading an element, damage 
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accumulates within the material. Miner and Palmgren have postulated that failure 
resulting from multistage loading can be defined as: 

n - 

(1) Sr s 1 Miner's Rule 

N i 

Miner's Rule assesses the "damage" produced by a stress a. for n • cycles as n^/N^ and 
assumes the individual damages are additive and independent of sequence. [36] 

If Miner's Rule is assumed to correctly predict the cumulative damage for the 
8MT, the following equation can be stated: 

(2) t i /T 1 + t 2 /T 2 + t 3 /T 3 + ... + t n /T n = 1 

where t = time at load F„ 
n n 

T = Total time to failure at given stress level (S-N curve) 
n n 


Since T = L , where L = life to failure predicted by a Wei bu 1 distribution. 
Equation (2) n can b^ rewritteR as 



or (4) 

Harris [1] defines a relationship based on the dynamic loading capacity as 



where C = dynamic loading capacity 
P = 3 for rolling bearings. 




allows the life to be related to the applied axial reaction load and to failure at a 
given stress (force) level. 


Defining Equation 6 in general terms gives 


(7) 


'n + 1 




Substituting Equation (7) into (3) and assuming 




'77W 

( 8 ) 

A similar equation can be written for L 2 such that 


3-stage loading gives 
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or 



Equation (10) may be expanded to a general form: 



Recall that a summation can be represented by the integration of a function so 

that 



Therefore , 



the integral term represents the area under a (Fi) p vs. Time curve. If a plot of 

( Fl ) VS ,, tim , e J s - siven ( See Fi 9 ure 5 -l-Z), the life at a particular load can be 
graphically determined by calculating the area under the curve using the prescribed 
boundary conditions. 



FIGURE 5.1.2 LOAD VS TIME 
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5 . 1.3 Application of the Double Linear Damage Rule to Bear ing Life Predictions 

Althouqh Miner's Rule has been used to analyze cumulative damage, it does not 
account for the effects of load sequencing. Miner's Rule assumes that when the sum- 
mation of cycle ratios (n/N) equals one, failure occurs. How ®ver, the cycle ratios 
can sum to be less than or greater than one when high-low or low-high load sequenc 
ing, respectively, is followed. 

In the BMT, bearings will experience various load sequences and it is necessary 
to account for these associated effects to correctly predict bearing life. The 

application of the Double Linear Damage Rule (DBDR) to bearing life pred-ictions is 

being evaluated. Most of the work by Halford at the Lewis Research Center has been 
with 0.25 inch diameter specimens. Discussions have been held . Dr ; 
concerning the feasibility of applying the DLDR to bearings, and indications are that 

it can be applied. 

5 . 1.4 As sessment of 57 mm Turbine Bearing Life as a Functions of Ra dial and Axia l_ 
load - 

An analysis was conducted to evaluate the 99% (B.) survival life of the 57mm LOX 

pump turbine end bearing as a function of axial and radial load. The analysis 

consisted of using the B in (90% survival life) predicted by the Shaberth I bearing 

code and 
life. The 

was used to account for the . - — . - . 

cryogenic coolant. It is recognized that the bearing life formulas used in the 
computer codes are for conventionally lubricated bearings and extrapolation to the 
cryogenic environments experienced by the turbopump bearing produce unknown errors. 
However, because of the unavailability of life data for bearings operating in cryo- 
genics, these are the best estimates possible and should at least provide a measure 
of the sensitivity of bearing life to the load factors. 

In general, a Weibull distribution adequately correlates the life dispersion of 
conventionally operated rolling bearings experiencing fatigue type failures in we 
range of 10% to 50% failures. Designers of aerospace hardware are interested in 
reliabilities of 99% or greater. Consequently, the deviation of the available test 
data from the Weibull distribution is of particular interest in the region of high 
survival probabilities. The following is a synopsis of the method used to estimate 
the life information shown in Figure 5.1.3. 


i of using the b in survival me; 

modifying the results, as described in previous sections, to estimate the B. 
he lubrication life reduction factor used for the Shaberth IJ bearing code 
to account for the reduction in life due to marginal lubrication in the 




Assuming a Weibull distribution, the bearing life can be related to a given 
probability of survival as follows: 


1 ) 


B = B 


10 


fefJ 


1/e 


B bearing life 

B. n bearing life for 90% survival 

S 1U Probability of survival 

e dispersion parameter (10/9 for ball bearings) 
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Equation 1 does not adequately represent the test experience for high probabil- 
ity of survival (S>S0%). As an example, for 100% survival Equation 1 predicts zero 


life. A technique has been developed in 
between the Wei bull distribution and test 
developed in the reference. 

1/S + Ye 


Reference 44 to correct for the deviation 
experience. The following expression was 


10 


[ L n 1/S + Ye 
L n 1/.9 _ 


1/e 


2) B = B 

1/ . 9 

Ye is the standardized excess life as defined below: 


Probability of 
Survival % 


Standardized 
Theoretical Life 


Standardized Excess 
Life 


Y <_ .001 
.001 < Y < .05 
.05 < Y < .9 


Y + Ye = .004 

Ln ( Y + Ye) = 
0.69 Ln (.328Y) 
Ye = 0.013 


S 99.9 
99.9 ^ S _> 95 

95 i s L 40 

3) Ve = fe) eL » 2 

Equation 3 can be expressed in terms of survival probability as follows: 


4) Y = 


= / L n 1/S 
\n l/.l 


L n 2 ’ L n 


Therefore, for 99% survival 


Y s L n .99 " * 01 

The standardized excess life can be determined from the appropriate expression 
in Table 1. 


w 


6) Ye = [. 328Y j ' 69 - Y = [(.328) (.01)] * 69 - .01 = .0093 

The relation between and B 1Q life can be determined from equation 2. 
d _ n f L n 1/.99 + . 0093 ~ 1 * 9 . 

1 «L w'g J =0 - 22B io 

Without the . excess _ life term (Ye), the B. life estimate would be (.12 B. 0 ). 
Therefore for this specific case, the B. life is approximately twice that predicted 
by the Weibull distribution. The other factor considered in this analysis is the 
reduction in life due to marginal lubrication. Lubrication reduction factors used in 
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the Shaberth II bearing code are as low as 0.21. Applying this value to Equation 7 
gives: 

8) Bj = .046 B 1q 

The bearing life data shown in Figure 5.1.3 were generated by adjusting the B in 
life predicted by Shaberth I, as indicated by Equation 8. As shown, the bearing lire 
at the higher radial loads is relatively short. Due to the uncertainties in the 
extrapolation of life data from conventional to cryogenic environments, the results 
shown should be used with caution, especially in the regions of high radial loading 
(> 2000 lbs). Other factors such as contact temperature and dry lubrication perfor- 
mance at elevated temperature have not been considered. These factors will be more 
critical at the higher loads due to increased heat generated in the contacts. 

The following is a suggested technique for relating bearing life to contact 
temperature. This method relies on the assumption that the dominant temperature 
effect is the loss of hardness in the contacts. This assumption may be adequate up 
to the temperature level at which the solid lubricant begins to deteriorate. 

Several investigations (References 45 and 46) have attempted to relate bearing 
fatigue life to surface hardness. As pointed out in Reference 45, hardness appears 
to be a dominant factor determining life, all other factors being equal. Therefore, 
if data were available to relate bearing fatigue life to hardness, and hardness to 
temperature, the effect of bearing temperature could be factored into the life 
equation. This proposed procedure is shown graphically in F ; „.ure 5.1.4. This is an 
example only and does not reflect the performance of 440C material. As discussed 
later, the necessary data to construct such a chart for 440C has not been found in 
the literature. 


The reduction in bearing capacity and B.^life can be obtained from bearing 
fatigue testing at specific hardness values. ^The relationship between temperature 
and hardness is needed to factor temperature into the life equation. As an example 
(See Figure 5.1.4), if a bearing contact is operating at 400°F, the dynamic capacity 
would be reduced by a factor of 0.81 and the life would be reduced by a factor of 
0.56. This assumes the bearing life is related to the cube of the load. Review of 
the literature indicates a deficiency in the data necessary to construct curves for 
440C such as those shown in Figure 5.1.4. Fatigue data relating the life of 440C to 
hardness, and hardness to temperature, would be of considerable interest. The 
availability of this information would allow the temperature estimates from the 
bearing thermal model to be used to gain further understanding of the thermal effects 
on bearing life. 

5.1.5 Investigation of Ball Spin-to-Roll Ratio 

An important factor in bearing life is the amount of slip occurring between the 
ball and races. The Lundberg - Palmgren method (Reference 1, page 411) for evalu- 
ating rolling bearing fatigue life is based on the assumption that bearings are lu- 
bricated and operated "properly". Proper operation means proper alignment, adequate 
loading to preclude skidding, free of contaminants such as dirt, etc. When the 
bearing is operated properly, sliding occurs in the contact area due to deformation 
of the surfaces. Since the Lundberg- D almgren fatigue model is based on empirical 
data, this type of sliding is assumed to be accounted for in the fatigue life 
equations. Sliding which occurs because of skidding is not accounted for. According 


78 



spectra research systems 


original PAQE 

OF POOR QUALITY 


to Reference 1, the following expression can be used to reduce the dynamic capacity 
of the point contact due to ball spin. 


Q 


c 



_u Jfg. 


b 

a U 


roll 


* 

% 


Presumably 0 must be determined experimentally, and 


b = Semi minor axis of contact ellipse 

a = Semi major axis of contact ellipse 

W = Ball spin normal to contact surface 

W s ,, = Ball roll relative to contact surface, 
roll 

As a start in investigating this effect, the ball spin to roll ratio for the 57 
mm bearing was calculated as a function of AT across the bearing. Since the SHABERTH 
computer program does not provide this calculation directly, additional computations 
were performed. The necessary ball rotational information generated by SHABERTH is 
as follows: 


W = x component of ball speed 

W x = y component of ball speed 

\r = z component of ball speed 

W z = Resultant ball speed 

r 

Using this information, the ball spin (W ) and roll (W r ) can be determined by 
constructing a vector diagram as shown in Figure 5.1.5. From Figure 5.1.5, the 
following expressions can be obtained: 

Ball Spin: 

W si = Sin ct i +w x sin “i" w y Cos a i 

W - W « -W Sin a +W Sin a -W Cos a 
so so o ox o y o 


W ri ” (w s' w 0 > Cos Cos a.+W Sin a, 

u so IX i y 1 

w - W Cos a-W Cos a -W Sin a 
ro o ox o y o 

Where: W^, W so = Ball spin normal to inner and outer races respectively. 

W s = Shaft speed 

W = Ball orbital speed 

o 

a., a Q = Inner and outer contact angles 

W W - Ball roll speed relative to inner and outer races 

ri ro 
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The spin-to-roll ratios can therefore be expressed as follows: 


Inner Race: 


w 

w 


si 


ri 

Outer Race: 

w 


so 


V 


ro 


<W g -W o ) sin a in a^-W y cos 

(W-W) cos a.+W cos a.+W sin a. 
so i x i y i 

-W sin a +W sin a -W cos a 
o o x o y o 

W cos a -W cos a -W sin a 
o ox o y o 


The inner and outer race ball spin-to-roll ratios are shown in Figure 5.1.6 as a 
function of bearing temperature difference. As shown, the ball spin is reduced as 
the temperature difference is increased. This is the result of increased normal 
loads and friction forces and the reduced differences in inner and outer contact 
angles. It is evident that for a given set of conditions; load, speed, friction 
coefficients, etc., there is an optimum bearing configuration, contact angle, inner 
and outer race curvature, etc., that will minimize the combined effects of ball spin 
and heat generation and maximize bearing life. These analyses have been conducted 
using a clearance stack producing a minimum contact angle. 


5,2 Shaft Deflection and Bearing Reaction Analyses 

5.2.1 IPX Tester Shaft Deflection as a Function of Bearing Eccentricities and Radial 
Loads 

The results of the deflection analysis to support the redesign of the LOX tester 
are shown in Figure 5.2.1. These data were generated from a SHABERTH model of the 
LOX tester shaft bearing system shown in Figure 5.2.2. The deflection curve shown is 
the sum of deflections due to bearing eccentricities, outer race to carrier 
clearances, and bearing and shaft deflections due to radial loading. As shown, the 
"no load' 1 (deflection at zero radial load) is dependent on the assumed component 
operating temperature. Over the temperature range investigated, this is a small 
effect. All bearings except No. 4 was assumed to be axially perloaded to 1,000 lbs. 
Bearing No. 4 preload was 500 lbs. Five hundred pounds is the minimum preload 
expected with an axial load large enough to bottom the preload spring of bearing No. 
3. Since the bearing radial stiffness decreases with axial preload, this is a worst 
case conditions for shaft radial deflection. 


The individual bearing reactions as a function of radial load are shown in 
Figure 5.2.3. As shown, the applied load is not distributed equally to the bear- 
ings. The principal reasons for this are that the bearings are not located equi- 
distant from the point of load application and the axial preload is not the same for 
all bearings. These differences will change as the shaft is axially loaded and 
deflects in the axial direction. For these reasons, bearing radial reactions should 
be measured at each bearing to provide known loading conditions for correlation of 
loading effects on bearing life. 

5.2.2 Bearing Reactions and Deflections for the LOX Bearing Tester 

This analysis was done for the latest preliminary design of the LOX Bearing and 
Materials Tester. The model used is shown in Figure 5.2.4. The results of the 
deflection analysis are shown in Figure 5.2.5. The deflection curve shown is the sum 
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of deflection due to bearing eccentricities, outer race and carrier clearances, and 
bearing and shaft deflections due to radial loading only. All bearings were assumed 
to be axially preloaded to 1,000 lbs. As shown, the maximum shaft deflection for the 
4,000 lbs radial load was found to be approximately 2.5 mils. Figure 5.2.6 shows the 
deflections for each bearing on the shaft as a function of radial load. 

Figure 5.2.7 shows the individual bearing radial reactions as a function of 
applied radial load with no axial loading. An axial load will change the distribu- 
tion of these radial reactions and further analysis would be required to estimate the 
radial reactions under combined radial and axial loading. 

5.2.3 Bearing Axial Reactions Resulting from Carrier Axial Loads 

The design of the Bearing and Materials Tester requires that bearing axial loads 
be applied by loading the bearing carrier rather than the shaft. The bearing carrier 
can be loaded by activating the axial load cells. Another source of loads on the 
carrier is the fluid pressure difference across the carrier. The fluid pressure over 
the carrier faces is difficult to estimate because of the uncertainty in the fluid 
rotation effects on local static pressure; therefore, the resulting pressure load is 
uncertain. Figure 5.2.8 provides the relationship of bearing carrier loads to 
bearing reactions. Since it is difficult, as explained in Section 5.10, to separate 
the load cell and fluid AP loads, the carrier loads in Figure 5.2.8 should be 
considered the sum of these loads. As shown, the preload spring constant uncertainty 
affects the initial preload about 20% and has a smaller effect as the carrier load is 
increased. Also of interest is the reduction in preload of bearings 1 and 4. 
Although the preload for these bearings is reduced to about 400 pounds, it has been 
shown in a previous analysis (see Section 5.8) that this load is sufficient to 
prevent ball skid due to fluid drag. 

5.3 Bearing Stiffness Analyses 

5.3.1 Stiffness Characteristics of the Turbine End LOX Pump Bearing 

An analysis was conducted to estimate the radial stiffness of the LOX pump 
turbine end bearing as a function of axial preload and shaft speed. The bearing 
radial stiffness is defined as the ratio of radial load change to radial deflection 
change (dF/d6). 

Approximately 21 computer runs were made using the SHABERTH bearing analysis 
program. These runs were made to determine the relationship of load versus de- 
flection as affected by axial load and shaft speed. The bearing stiffness was then 
determined as the slope (dF/dfi) of the force versus deflection curves. The radial 
load was varied from 100 to 400 lbs. At higher loads these curves can become nonlin- 
ear, i.e., the deflection can decrease as the radial load increases. As shown in 
Figure 5.3.1, the relationship is linear for axial loads greater than 1136 lbs. For 
the axial load case of 607 lbs, the curve is linear to about 200 lbs radial load. 
Since the computer solution failed to converge at higher radial loads, it is not 
recommended to extrapolate the stiffness versus axial reaction curve shown in Figure 

5.3.2 to axial loads below about 600 lbs and radial loads greater than 200 lbs when 
the axial reaction is less than approximately 1000 lbs. 


As expected, the radial stiffness of the bearing increased with increasing axial 
reaction or preload as shown in Figure 5.3.2. As previously discussed, the curve 
shown in Figure 5.3.2 represents the linear portion of the radial reaction versus 
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FIGURE 5.2.7 APPLIED LOAD VS RADIAL REACTION 
(LO* TESTER) 


30,000 RPH 

T * -IB0“C (NO AT ACROSS BEARINGS) 
NO AXIAL LOAD 

ALL BEARING PRELOADS = 1000 LB 
TESTER DIMENSIONS 
SCALED FROM DRAWING 
A30A852Q0- 1 


BEARING RADIAL REACTION. LB 


57nm BEARING 
30,000 RMP 

NO AT ACROSS BEARING 
-180°C TEMP 


BEARINGS 2 A 3 

SPRING CONSTANT « 55 Lbs /mill 

SPRING CONSTANT » 65 Lbs/mill 


BEARINGS 1 b 4 


CARRIER AXIAL LOAD X I0 5 Lbs 


FIGURE 5.2.8 BEARING REACTION Vs CARRIER AXIAL LOAD 
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deflection bearing characteristics shown in Figure 5.3.1. A similar analysis was 
done to determine the effects of shaft speed on bearing stiffness. Figure 5.3.3 
provides the bearing radial deflection as affected by radial reaction and shaft 
speed. Shown in Figure 5.3.4 is the bearing radial stiffness as a function of shaft 
speed. As shown, the bearing becomes less stiff in the radial direction as the shaft 
speed increases. This is primarily due to the increase in centrifugal force effects 
on the outer ring. 

5.3.2 Bearing Radial Stiffness as a Function of Ball Wear 

The SHABERTH I computer code was used to evaluate bearing radial stiffness as a 
function of ball wear. Although the highest level solution would not converge, a 
lower level solution that still considered friction was obtained. Judging by the few 
third level computer runs that did converge, the calculated radial stiffness of the 
bearing did not differ significantly between solution levels. 

Four radial loads {100, 200, 300, and 500 Tbs) were applied to bearings at five 
different stages of ball wear. The results of these runs are shown in Figure 5.3.5. 
This analysis suggests that even a decrease in ball diameter of 0.0008 inches is 
enough to reduce the radial stiffness by half. 

5.3.3 Investigation of the Radial Stiffness of the BMT 57 mm Bearings 

An analysis has been conducted to evaluate the effects on bearing operating 
characteristics (stresses, radial stiffness) caused by impressed radial loads. The 
purpose of this investigation is to identify these characteristics tor input to test 
planning to determine the radial dynamic response of the bearing shaft system. The 
radial load range investigated was 500 to 6000 lbs. with axial loads of approximately 
500 to 2800 lbs. Assuming an impressed axial load of approximately 2500 lbs due to 
the pressure differential across the bearing carriers, the two inboard bearings will 
be axially loaded to that value and the two outboard bearings will be unloaded. 
Preload is 1000 lbs to approximately 500 to 800 lbs. Since axial loads tend to 
stiffen an angular contact ball bearing in the radial direction, it was expected that 
the unloaded bearing would be the most distressed at high radial loads. 

Shown in Figure 5.3.6 are the bearing deflections versus radial reactions. As 
expected, the bearing with the larger axial load exhibits the larger stiffness. 
There is also a stiffness change as the load increases. This is due to the increased 
radial deflections unloading balls, i.e., the total load is supported by fewer 
rolling elements. It is interesting to note that the unloaded bearing has three 
apparent stiffness values. The first change in slope (stiffness) occurs for the same 
reason as described above. The second change in stiffness is believed to be caused 
by the remaining load carrying balls being forced down on the inner race shoulder. 
Since at this location the inner contact angle is approaching zero, the only radial 
deflection possible is caused by compressing the balls (i.e., all ball movement is 
prevented). Shown in Figure 5.3.7 is the inner contact angle as a function of radial 
loads for the unloaded bearing. The bearing is reacting the radial load at an 
azimuth angle of zero. Notice that at this position, the inner contact angle dras- 
tically decreases as the load increases. As this happens, the ball is forced down to 
the shoulder of the inner race. Also observe, at an azimuth angle of 166°, the inner 
contact angle goes to zero for radial loads above about 1000 lbs. This indicates 
that the balls at this location have left the inner race and, due to centrifugal 
force, are following the outer race. At the higher radial loads only three balls are 
supporting the load. At the higher axial loads, this does not occur, and all balls 
are in contact with the inner race. 
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growth. The following conditions were assumed in the thermal analysis: 

Inlet fluid temperature -305°F 

IP?* . , 6-3 lbs/ sec per bearing set 

Axial load 2500 lbs, 

ahoMt tS ? 7 or dl -^ 6C L an average temperature difference between carrier and housing of 
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As clearances will exist during the early part of the test, analyses were done 
to investigate the effects of misalignment on bearing operation characteristics. 
Angular misalignments were determined as a function of carrier to housing clearances 
and the effect of these misalignments on bearing radial stiffness, ball speeds, 
contact angles, contact stress, and contact heat generation was determined. These 
effects were investigated for axial loads from 500 Ibs/brg to 6000 lbs/brg, and 
radial loads of 100 to 500 lbs. 

Figure 5.4.1 shows the potential for the carrier to tilt in the housing. The 
tilt angle 0 is a function of the carrier diametrical clearance (C). From the 
geometry of the figure the following relations can be written: 

1) M = W Cos e - D Sin 9 

2) M = V D 2 + W 2 - (0 + C) 2 

3} D + C - D Cos e + W Sin e 

Eliminating M and Cos e from the above: 

4) 0 = Sin" 1 f W (D + C) - dVd 2 + W 2 - (D + C) 2 * 

l o 2 + m 2 J 

Equation 4 provides the tilt or misalignment angle as a function of carrier dia- 
metrical clearance. Figure 5.4.2 provides misalignment angle as a function of 
diametrical clearance. As indicated the maximum expected clearance is 1.3 mils 
(carrier and housing at same temperature) giving a maximum misalignment angle of 2.25 
minutes (.0375°). 

5.4.1 Bearing Radial Stiffness Variations 


The bearing radial stiffness as a function of azimuth angle is shown in Figure 
5.4.3. In this figure, the zero azimuth angle corresponds to the "pinch" point 
caused by the misalignment of the outer race [the results from the "SHABERTH" bearing 
analysis program. Figure 5,4.4, are presented with the zero azimuth angle 180° out.of 
phase from the "pinch" point]. The variation in stiffness shown is conservative 
because the reduction of the bearing misalignment angle as a function of axial load 
was not considered in the stiffness determination. These variations should not 
significantly affect synchronous radial motion of the shaft due to small rotor 
unbalance which was the main concern of possible large stiffness variations due to 
misalignment. 

5.4.2 Bearing Oper ati on Characteristics 

The investigation of carrier misalignment on contact stresses, contact angles, 
ball speeds, and heat generated at the contacts required numerous "SHABERTH" computer 
runs with varying input parameters. The conditions shown in Figure 5.4.4 were 
investigated. 

In reviewing the data, the most significant variations due to outer race mis- 
alignment is the ball orbital speed. This occurs for the 500 lb axial load condi- 
tion. This is representative of tester bearings 1 and 4 with the normally accepted 
axial PA load of 2500 lbs. Comparisons of contact stresses, ball speed, and contact 
angles for the misaligned and normal conditions do not show significant variations in 
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FIGURE 5.<S.1 CARRIER TILT IN HOUSING 
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FIGURE 5. A. 2 OUTER RACE MISALIGNMENT AS A FUNCTION OF BEARING 
CARRIER DIAMETRICAL CLEARANCE 
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FIGURE 5.4.3 RADIAL STIFFNESS VERSUS AZIMUTH ANGLE 
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FIGURE 5.4.4. CONDITIONS INVESTIGATED FOR BEARING MISALIGNMENT EFFECTS 
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these parameters for the higher axial load cases. There is a maximum difference in 
ball orbiting speed of about 20 rad/sec for the 500 lb axial load condition. It is 
estimated that this causes + 8.5 mil movement of the ball relative to the ball pocket 
center. This should not cause a problem since there is 25 mils diametrical clearance 
between ball and ball pocket. 

5.4.3 Effect of Axial Load on Bearing Misalignment 

The effect of axial load on operating misalignment is shown in Figure 5.4.5. As 
the axial load is increased, the inner race is forced to rotate thus decreasing the 
angle between the inner and outer race. The degree of rotation is obviously depen- 
dent not only on bearing characteristics but also the stiffness of the shaft since 
the shaft exerts a clockwise movement on the bearing inner race as the race is forced 
to rotate counter-clockwise due to the outer race misalignment. Shaft and bearing 
characteristics were included in the bearing/shaft model. 

In summary, the analyses indicate no tester operational problems due to outer 
race misalignment with a carrier to housing clearance of 1.3 mils. Although the 
thermal analysis indicates a slight carrier to housing interference for steady state 
conditions, verification of this condition is necessary during early tester operation 
and should not be detrimental during the initial test phase in LN„ when loads are 
relatively light. For higher loads and increased heat generation, Turther analysis 
is required to evaluate the adequacy of carrier to housing clearance. These analyses 
should use the thermal model "calibrated" by thermal data taken in the early BMT 
runs. 


5.4.4 Potential Effects of Outer Race Misalignment on Shaft Assembly Balancing 

An apparent shaft assembly unbalance was observed during balancing of the BMT 
rotating assembly. The cause was determined to be bearing carrier face run out. The 
carrier run out can in effect produce an angular misalignment in the bearing outer 
race. The purpose of this investigation is to develop an explanation as to how the 
outer race tilt could produce a synchronous load effect, which would appear to be a 
grossly out of balance shaft assembly. 

As a start, the geometry change in the bearing was evaluated as the outer race 
was displaced. Figure 5.4.6 illustrates the geometric representation. As the outer 
race is rotated counter clockwise, the right ball moves upward and to the right and 
the left ball moves down and to the right. The center of the pitch circle is dis- 
placed in the positive y direction. The center of the circle formed by the outer 
race curvative center is displaced in the y direction (k) and in the x direction (h). 
Normally this is a static condition and the ball pitch circle would be displaced from 
the shaft center line by a constant value and there should be no shaft dynamic 
effects. However, the contact angle changes as a function of the azimuth angle. The 
contact angle is increased in the ^ = G° position and decreases when = 180°. At \p 
= 90 the contact angle is not changed. Assume an out of balance shaft assembly 
producing a synchronous load in the radial direction. Since the contact angle 
changes with azimuth location, the bearing radial stiffness also changes in the 
azimuth direction. Therefore, with a variable stiffness, the shaft can move radially 
in a dynamic manner synchronous with shaft speed. Consequently, a small unbalance 
can be magnified into an apparent large unbalance by outer race misalignment. The 
geometry shown in Figure 5.4.6 is obviously for an unloaded, static bearing. Loads 
and speed will change the configuration. However by superimposing the outer race 
misalignment effects on the configuration assumed by the bearing at speed and under 
load, the above explanation should still be valid. 
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To investigate the sensitivity of the contact angle to outer race misalignment 
the following analysis was made. Assuming no deformation in the contact areas and 
very small axial displacement, h^o, the contact angle can be related to the misalign- 
ment angle as follows: 


SINB ‘ = R o 6 COSt + $INB„ 

R 0 

A is the distance between curvature centers. 

The above equation is presented graphically in Figure 5.4.7. A one minute angle 
of misalignment introduces about 4.5% peak to peak variation in contact angle. The 
degree this affects radial stiffness and shaft assembly dynamics is yet to be de- 
termined. 

An analysis to determine the stiffness variation is in work. "SHABERTH I" was 
coded to evaluate the stiffness variation with misalignment, but failed to run. 
"SHABERTH 11" was subsequently coded to include outer race misalignment and it ran. 
These results will be bearing radial stiffness values that can be input to a bearing 
dynamic model for evaluating rotor/shaft dynamics. 

5 • s Bearing W ear Track, Ball Pass Count, and Ball Spin Speed Analysis 

A 57 mm bearing was modeled on the SHABERTH computer program with maximum clear- 
ances and an arbitrary temperature distribution radially across the bearing. The 
purpose was to determine if the observed bearing wear tracks could be matched within 
the known loads applied to the bearing tester. As shown in Figure 5.5.1, a fairly 
good match can be obtained with an axial load of 6500 pounds. For comparison, the 
calculated tracks for a uniform bearing temperature and nominal clearances are shown. 
The temperature gradient used was an assumed one which was selected to investigate 
the sensitivity of contact angle to increased diametrical clearances due to tempera- 
ture differences. This information shows that the bearing tracks can be influenced 
by temperature and clearance as well as applied loads. 

Cage speed data generated from SHABERTH was used to generate inner and outer 
race ball passes per minute. This data was compared with the data measured during 
tester operation. The comparisons, shown in Figure 5.5.2, indicate good agreement 
with the measured and calculated values. Also shown is the ball spin speed compared 
with the data derived from the tests. As illustrated, the above parameters are 
plotted as a function of the bearing axial reaction load. 

5.6 Effects of Ball Drag on Operating Characteristics 

The SHABERTH computer program has the capability of evaluating the effects of 
ball drag. Unfortunately , the treatment of the hydrodynamic effects have been 
programmed considering the properties of lubricants operating in a nominal tem- 
perature range of -60 to 600°F. For example, the temperature correction for the 
lubricant density is given as: 


p (T) = p (60°F) - G(T - 60°F) 

where p (60) and G can be input or a lubricant can be selected that has its properties 
already programmed. Obviously, 0~ does not exist as a liquid at 60°F and the above 
equation does not account for a phase change from gas to liquid. Therefore, a ficti- 
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tious value for p(60) can be used or a lubricant selected that has a density range 
covering LOX, when the appropriate temperature is used in the above equation. 

Another problem associated with LOX as a fluid is caused by the way the program 
accounts for viscosity variations with temperature. Typical lubrication viscosity 
data have been correlated with the following expression: 

log 1Q log 10 (V + 0.6) = A - B log 10 (T + 460) 

where A and B are constants determined by substituting the known values of viscosity 
(V) at T = 100°F and T = 210°F and solving the two resulting equations for A and B. 
For a low viscosity fluid such as LOX, the term (V+.6) is less than one thus re- 
quiring taking the log of a negative number. 

To work around the above problems, SHABERTH was run at LOX temperature using a 
dry lubricant friction coefficient. The purpose was to determine the operating 
clearances at LOX temperature. These clearances will be used as input for a higher 
temperature run (temperature selected to obtain a match of lubricant and LOX density) 
with a fluid to determine the drag force on the balls. 

T he ball drag force is determined by: 

Fw = p gg — (dmWo) 2 

A v = ball frontal area 

Cv = drag force; determined as a function of Reynolds Number 
dm = bearing pitch diameter 
Wo = ball orbital velocity (cage speed) 
p = density of fluid in the bearing cavity 

The density of the fluid in the bearing cavity is determined by: 

P = XCAVpo 

p is the density of the liquid and XCAV is an input variable defining the percent 
l?quid in the bearing cavity. Obviously, the ball drag force can be varied by 
varying the value of XCAV. The technique for using SHABERTH for estimating balT drag 
effects is as follows: -*• 

1) Run SHABERTH at LOX temperatures to establish operating clearances. 

2) Run SHABERTH at appropriate temperature to provide a lubricant density 
equivalent to LCX density. 

3) Use the values of viscosity to calculate a Reynolds Number to determine the 
drag coefficient Cv, and calculate the drag force used by the program. Any 
value of drag force can be investigated by varying XCAV. 

To date, successful runs have been made using this technique. The results are 
provided in Section 5.7. 
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5 . 7 Estimates of Ball Drag Forces, Torques, and Heat Generation 

Proper manipulation of input data allows the SHABERTH bearing program to calcu- 
late the heat generated due to viscous work and drag on the orbiting balls in a low 
viscosity fluid such as LOX. This technique was explained in Section 3.4.6. Esti- 
mates of ball drag, torque, and heat generation can be made as a function of the 
percent liquid assumed to be in the bearing cavity. Shown in Figure 5.7.1 are the 
preliminary results of an analysis on the 57 mm bearing. As noted, the data are for 
a single bearing assumed to be operating in a stagnant fluid. The percent fluid in 
the cavity can be interpreted either as a low density fluid or the partial filling of 
the cavity with liquid. 

Although bearing characteristics such as contact angle were not affected by ball 
drag, the fluid drag on the balls does cause ball slippage or skidding in a lightly 
loaded bearing. This effect is shown in Figure 5.7 . 2 which gives the ratio of 
cage-to- shaft speed as a function of bearing reaction, percent fluid in the cavity, 
and friction coefficient. Although this information should be considered preliminary 
due to convergence problems with the machine solution, the trends are believed to 
represent the physical characteristics. This problem should be overcome by going to 
a higher level of solution in the SHABERTH program. 

As shown, the effect of increased flooding of the bearing cavity is to slightly 
reduce the cage rotational speed at bearing reaction loads above approximately 150 
pounds. Below 150 pounds there is a significant reduction in cage speed indicating 
bearing slip on the inner race. A reduction in friction coefficient increases the 
ball slip which further reduces the cage speed as shown. This indicated that a com- 
bination of low axial load and high fluid drag can cause significant slip between the 
rolling elements and the inner race. The consequence of this is increased heat gen- 
erated, wear, and shorter bearing life. The fluid density used in these analyses is 
representative of LOX. 

The effects on bearing frictional and viscous heat generation are shown in Fig- 
ures 5.7.3 and 5.7.4. Figure 5.7.3 provides data for a 50% flooded cavity and Figure 
5.7.4 is data representing a 100% flooded cavity. As shown, the fractional heat 
generated increases, after passing through a minimum point, as the bearing reaction 
load is reduced. The effect is more severe for the 100% flooded cavity. _ These 
effects could be considerably more severe than indicated for some combinations of 
axial and radial loading. For instance, if a bearing is lightly loaded in the axial 
direction with significant slip between rolling elements and inner race, a sudden 
application of radial load will cause the friction forces between the ball and race 
to increase. Since the ball/cage mass cannot be instantaneously accelerated to the 
cage equilibrium speed, considerable damage can be done to the balls and races due to 
slippage under high loading. This input could explain the abnormal wear observed in 
the LN 0 tests for bearings 1 and 4 of the BMT. Efforts are under way to extend these 
analyses to include other friction factors and modify the fluid density to represent 

ln 2 . 

5.8 Effects of Bearing Preload on Cage Slip and Frictional Heat Generation 

An analysis has been conducted to evaluate the effects of axial load on bearing 
cage speed and heat generation rates due to rolling element slip with respect to the 
inner race. The effects of viscous drag, due to the rolling elements rotating in the 
LN 0 coolant, was accounted for in modeling the 57 mm LOX turbopump bearing for 
analysis using the SHABERTH Bearing Analysis Computer Program. 
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The objective was to determine the axial load that would allow the viscous drag 
force to overcome the contact friction forces to the extent that significant slippage 
could occur at the rolling element/face interfaces causing appreciable increases in 
inner race heat generation. Since the friction force is directly proportional to the 
coefficient of friction, which is not conclusively known, the friction coefficient 
was varied over a range from 0.08 to 0.22. The results of the investigation are 
shown in Figures 5.8.1 and 5.8.2. The data shown in Figure 5.8.1 for a friction 
coefficient of .08 indicate that, below a specific value of axial load, the heat 
generated at the inner race decreases to a minimum and then increases to a maximum. 
The same characteristic should occur for the other friction coefficients. 

Inability of the computer program to converge to a solution at the lower levels 
prevented completing these curves in the low load region. It should also be observed 
for the friction coefficient of .08 that the heat generated by cage drag begins to 
decrease as the heat generated at the inner race begins to increase. Shown in Figure 
5.8.2 is the cage/shaft speed ratio as a function of axial reaction, friction factor, 
and fluid drag. The upper horizontal line represents the case for no fluid drag. 
For these conditions, the cage to shaft speed ratio is about 43%. As shown, a de- 
crease in friction force (friction coefficient) reduces the cage speed. Furthermore, 
a reduction in axial load reduces the cage speed dramatically at loads below about 
300 lbs. 

The results of this analysis provide a better understanding of why bearings 1 
and 4 of the BMT were degraded even though they were lightly loaded in the axial 
direction. Based on these results, a design requirement should be to provide at 
least 400 lbs of preload in the bearings at all operating conditions. 

5.9 Assessment of Bearing Wear by Measurement of Axial Deflection and Applied Load 

It is desirable to have the capability to evaluate bearing wear subsequent to 
each test run. In principal this can be done by measuring the shaft deflection as a 
function of applied axial load, and correlating this information with the analytical- 
ly developed curves shown in Figure 5.9.1. To assess the degree of wear, the infor- 
mation shown in Figure 5.9.1 was generated by developing the relationship between 
carrier axial load and bearing reaction and, by use of the Shaberth bearing code, 
developing the relationship between bearing reaction and shaft deflection. The data 
shown are for cryogenic temperatures since the deflection sensors are expected to be 
calibrated in this temperature range. Ambient temperature curves can be readily 
generated if required. 

5. 10 Bearing Operating Preloads 

An analysis has been done to estimate the change in bearing preload, from the 
assembled room temperature value, as the bearing/shaft assembly is cooled from room 
temperature to LN ? temperature and rotated to a speed of 30,000 RPM. At this point 
in the analysis, Tt is assumed that the shaft is not axially loaded; therefore, the 
bearing axial reactions are equal for all bearings and equal to the spring preload 
value at operating conditions. 

Considering only the thermally induced axial loads and deflection, the following 
expressions can be written for the axial deflection of the shafts and components. 


Heat Generated # Inner Race Watts 
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c - Components on the shaft 

Since di,,. = d£ these equations can be solved for the thermally induced force 
and the axial deflections due to the temperature changes. The axial deflection of 
interest is the net deflection between one-half the bearing separator, one-half the 
bearing race, the bearing spring, and the carrier tang. Using the data provided in 
Figure 5.10.1, a net deflection of 0.1551 mils was estimated. This deflection is in 
the direction to compress the spring. Since the spring constant is between 55 and 65 
Ibs/mil (depending on direction of spring travel), this movement produces an insig- 
nificant change in preload. 

An axial movement of the bearing can occur as a result of radial clearance 
changes caused by initial press fits, temperature changes, and shaft speeds. To 
assess these effects on bearing preload, bearing deflection curves at these different 
conditions were generated, as shown in Figure 5.10.2. Deflection curves are shown 
for three conditions; room temperature, cryogenic temperature, and cryogenic tempera- 
ture and speed. As shown, fitting the bearing to the shaft causes an axial movement 
of the races which requires about a thousand pounds to return the races to their 
relative unmounted positions. Since the room temperature preload is approximately 
one thousand pounds, this is the initial condition from which the shaft/bearing 
system changes as the components are cooled down and brought up to speed. The effect 
of the cool down is shown, as well as the speed effect. 

As indicated in Figure 5.10.2, the effect of cool down on bearing deflection is 
to release a small amount of preload, ^50 lbs. As shown in Figure 5.10.1, the 
thermal coefficient of expansion for the shaft is about 10% larger than the value for 
the inner race; therefore, the shaft shrinks from the inner race reducing the inter- 
ference fit. This, in turn, increases the internal clearance of the bearing, arflow- 
ing the outer race to move in a direction to reduce the spring preload. As the shaft 
reaches operational speed (30,000 RPM) , the centrifugal force causes the inner race 
to grow, reducing the bearing internal clearance, causing the outer race to move in 
the direction to increase the preload. As shown in Figure 5.10.2, these are compen- 
sating effects and the operating preload is very near the room temperature value. It 
is important to realize that these conditions have been estimated for uniform temper- 
ature conditions (i.e., all bearing elements, shaft, and carrier are assumed to be at 
the same temperature). Temperature gradients across the bearing and in the bearing 
carrier could change these results. In addition, the spring constant is the room 
temperature value. Further analyses are needed to evaluate these effects. 

Another factor affecting bearing preload is the preload spring hysteresis. This 
effect is shown in Figure 5.10.3. As shown, the spring has a different spring 
constant depending on the direction of travel. It is estimated that the preload 
value will vary about 20% depending on the direction of the spring action. 
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5.11 Interface Fit B etween the Inner Race and Shaft 

The bearing inner race is a force or interference fit on the tester shaft. This 
assembly is made at room temperature and the interference between the race and shaft 
changes as the tester is cooled to LNp temperature, rotated to test speed, and 
loaded. An analysis was done to estimate these changes at the various conditions 
from assembly to test. 

Assuming that the assembly temperature is 20°C and the cold operating tempera- 
ture is -180 C, the interference will change from 2.3 mils at room temperature to 
1.2074 mils at -180 C. Rotation of the bearing shaft assembly to 30,000 per minute 
will reduce the interference to 0.6178 mils. A preload of 1000 lbs. will increase 
the interference to about 0.75 mils. Increased axial loading will "tighten" the 
interference of bearings 2 and 3 and loosen bearings 1 and 4. The interference for 1 
and 4 will not be less than 0.6178 mils. 

The above estimates are for uniform temperature of bearing race and shaft. This 
is not the case when bearing friction heat warms the inner race to a higher tempera- 
ture than the shaft. The maximum temperature difference between the inner race and 
the shaft probably occurs during tester start when the heat generation, due to speed 
and load, is suddenly imposed on the chilled shaft and bearing. This transient 
condition has not been analyzed; however, a calculation was made to estimate the 
temperature difference between the inner race and shaft that would eliminate the 
interfergnce fit. Assuming the shaft remained at -180°C, this analysis indicated 
about 38 C temperature difference between bearing inner race and shaft will cause the 
interference fit to be lost (i.e., if the inner race temperature attains an average 
temperature of -142 C, while the shaft average temperature remains at -180°C, the 
interference between inner race and shaft will be zero at 30,000 RPM and 1000 lbs 
axial reaction at the bearing}. 

Since the thermal model predicts an average temperature of only 28°C for an 
axial load of 2500 lbs and 30,000 RPM, there should be no problem in losing interfer- 
ence fit at steady operating conditions. However, as stated previously, the start 
transient could cause a temperature difference large enough to eliminate the inter- 
ference fit for short periods of time. Although an increase in axial load tends to 
increase the interface pressure (between inner race and shaft), it also increases the 
heat generation rate. The steady conditions analyzed were for an axial load of 2500 
pounds. It is expected, but not confirmed by analysis, that the increase in load 
will outweigh the thermal effects and the interference fit will be maintained at the 
higher load conditions. If the maintenance of an interference fit is critical to 
tester operation, a transient analysis should be done and additional higher loads 
investigated to verify that the interference is not lost for any possible operatinq 
condition. 

The previous analyses show that the room temperature bearing preload is not 
significantly affected by cooling to LN„ temperature and rotation of the shaft to 
30,000 RPM. Temperature gradients across the bearing have been neglected and the 
temperature effect on bearing spring constant has not been included. Although the 
bearing internal clearances can be significantly affected, it is expected that 
temperature differences of less than about 150°F across the bearing will have little 
effect on the preload value. The temperature effect on spring constant will require 
further work. It should also be observed that the effect of speed is to increase the 
preload; therefore, at a shaft speed of 10,000 RPM, the preload will be slightly 
reduced (reduction will be less than ^65 lbs.) compared to the value at 30,000 RPM. 
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The analyses relating the carrier load to bearing reactions shows that the un- 
loaded bearings (1 and 4) will retain sufficient preload to prevent excess skidding 
and internal heat generation. It also shows that the uncertainty in bearing reaction 
due to preload spring hysteresis is approximately 20% maximum. 

The interference fit between the inner race and shaft has been evaluated to 
determine the effects of temperature and speed. Although an interference is main- 
tained at uniform operating temperatures and speed, temperature differences between 
the bearing inner race and shaft can be large enough to eliminate the interference. 
It has been estimated that interference fit will be lost at sp^ed if the inner race 
temperature is greater than the shaft temperature by about 38 C. This assumes the 
shaft to be operating at -180°C which is a representative temperature when operating 
in LN~. It presently appears that this can possibly occur only during start tran- 
sients when the shaft and bearing are chilled and a sudden load and speed is imposed 
which heats the bearing before the shaft temperature has time to respond. 

5.12 Fluid Film Thickness Analyses 

5.12.1 Estimation of Film Thickness for 57 mm Bearing Operating in LfL and LOX 

A cursory analysis was done to provide an "order of magnitude" estimate of the 
fluid film thickness for the 57mm bearing operating in LN^. Although the viscosity 
of fL is over two orders of magnitude less than typical lubricating oil, the high 
speea of the bearing could partially compensate and provide some measure of sepa- 
ration of the contact surfaces. This could be an aid to the solid film lubricant in 
providing contact lubrication. 

There are several major drawbacks in using current techniques to estimate fluid 
film thickness in N ? . These techniques have been developed for lubricating oils and 
such characteristics as compressibility, phase type, and extreme variations in 
temperature have not been considered. Furthermore, there is a lack of viscosity data 
for PL at the high pressures encountered in the contacts. Therefore, the error due 
to extrapolation to these pressures is not known. The assumption of incompressibil- 
ity should provide optimistic values of film thickness. The assumed variation of 
viscosity with pressure (exponential law) could also be optimistic. 

The results of the film thickness estimates are shown in Figure 5.12.1. As the 
temperature increases, the film thickness decreases rather suddenly and then levels 
off. The sudden decrease is characteristic of a liquid, and the leveling off is 
similar to the behavior of a gas whose viscosity increases with temperature. The 
estimate was not carried past 1000°R since viscosity for higher temperatures was not 
available. Figure 5.12.2 provides viscosity data at specified temperatures and pres- 
sures. Also shown are values for the pressure coefficient (x) and y'. The signifi- 
cance of these parameters will be discussed later. Comparing y' with the viscosity 
values at 400 psia shows fairly good agreement at the lower and nigher temperatures. 
Since the critical temperature of PL is 227°R, it is expected that the larger dis- 
crepancy at 300°R is due to rapidly ^changing properties near the critical point. 
properties were taken from Reference 35. 

As shown, the film thickness varies from about 1.8y to .03y inches. Since the 
surface finish of the balls and races are 0.5y and 6y inches respectively, the ratio 
of film thickness to aspherity height varies from 0.3 to 0.05. Furthermore, thermal 
analyses of the bearings indicate that the race way temperatures are considerably 
greater than 170°R. Consequently, the film thickness to roughness ratio would be 
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considerably less than 0.3. Bearing life adjustment factors, due to lubrication, 
does not adversely affect life predications for film thickness to roughness ratios 
less than about 4. SHABERTH I assumes dry friction for ratios less than 0.4. 

Following is an outline of the method used to provide the estimated film thick- 
ness data shown in Figure 5.12,1. 

Gurbin's formula for fluid film thickness is (see Reference 1): 


h° = 


1.95 R (GU)° ,7?7 
—0.091 


G = XE' 

Q z ■ 

E' = “2 
1-v 

(For Rolling Elements 
of Same Material) 


X= Pressure Coefficient of Viscosity 

(1 0 - Fluid Viscosity as Reference Pressure 

R = Equivalent Radius 

U = Entrainment Velocity 
= (uj + u 2 )/2 

E = Modulus of Elasticity 
v = Poisson's Ratio 


This expression was developed from the incompressible Reynolds equation and 
elastohydrodynamic considerations including effects of contact surface deformation. 
The fluid viscosity is assumed to vary exponentially with pressure, and temperature 
is assumed to be constant. As given above, the expression is for line contact and 
must be slightly modified to represent ball bearings (point contact}. 


The terms in the above expressions are developed and grouped as follows: 


Entrainment Velocity 



W g = Shaft Speed 
W c - Cage Speed 
o = Contact Angle 
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To determine the fluid entrainment velocity U refer to the above figure. 

Uj = 1/2 (dm - DCosct) (W -W ) 

U 2 = g W R Cosa 

For pure rolling, - U^, and 

U « 1/2 (H s - W c ) [^2] DCosa ; 7= ^ Cosa 

to express W in terms of W note that the cage velocity V is 1/2 the absolute shaft 
velocity V 1 , i.e. s c 


--V c - 1/2 dmW c = 1/2 V.i 

V. » W s (dm/2 - D/2 Cosa) = 1/2 W dm (1-y) 
1/2 dmW = 1/2 1/2 W dm (1-y) 

V J 


W c = 1/2 W $ { 1-y) 

Substitution of in the expression for U gives: 

u ■ i «s f-f 1 ] M C05 “ 

Therefore, the entrainment fluid velocity can be expressed in terms of bearing geome- 
try, contact angle, and shaft speed. 


Equivalent Radius (R) : 



_ R 1 R 2 _ 1/2D (1 -7) Cosa 

■ " rr-'T ) cosoi +7" 


For small values of a, Cos a ^1, and 


R = 1/2 D (1 - y) 


Load Per Unit Length : 

As stated earlier, the expression for Q 7 is for line contact and, for use with 
ball bearings, should be modified to represent point contact. The following repre- 
sents the two contact configurations. 
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The maximum Hertz Pressure for point contact and line contact are as follows: 

3 Q z 

CT mp = 2HF (Po-int Contact ) 

2 Q 2 

ff mf = "iTB" (Line Contact) 


Assuming equal Hertz pressures and loads for both contacts gives an equivalent lenqth 
s, f for the point contact. a 

o _ 4a 
E " T 

The dimensionless terms in Curbin' s equation can now be written as: 


G = 
0 = 



Cosa 


_ 3 Qz 

’ 7 aT"D{i-7) 



and Gurbin's equation can be expressed as: 

r -I 1 - 091 r 

ii° = 0.228 jl/2 D (1 -7 )| Tw s 

where W g = shaft speed in RPM 

A = pressure coefficient (in 2 /lb) 
p 0 ‘ = viscosity (Ibs-sec/in 2 ) 

E' = 33 x 10 ® psi for steel bodies 



, . This _ arrangement of variables illustrates the effects of speed, viscositi;. and 

oad on film thickness. The speed and viscosity parameters are raised to the same 
P™?r ar f influential. The lower power value for 

Aithntmh h t i ? e renders the fllm thickness rather insensitive to load 

inwS+ h 4 .l 9h , s . ha , ft spe ® d £ an offset low fluid viscosity, the viscosity for 

f?ym t thi t ckne e ss hl9h ° f th ® t6Ster iS St111 1nsufficient to Provide an 


the term 
changes. 
N~ is so 
adequate 


Estimation of Initial Viscosity and Pressure Coefficient (x): 

as follows* POnential ldW relating vis cosity to pressure used in Gurbin's equation is 


XP 

#*= e 
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implies that at P=0, p - p . Since the internal pressure of the bearinq 
p r e s s u re 1 $ 1 eve*! r0X7niate ^ 400 psia * a correction was made to p Q to account for this 


M= n Q e 


x(P - Po) = (mo e -> p °) e >p 


The parameter p' 
thickness expression. 0 
determined as follows: 


= tw 0 e AP0 ) is used as the ambient 
The pressure coefficient for each 


1 


viscosity in the film 
temperature level was 


(pp-pp >" 


Pj and P £ are stlown in Figure 5. 12.2. The estimated values 


The values of p., p„, r. auu r« 
of A and p Q are also given ifi this figure. 

h oa ^?n S l m j 1ar + * ana1 ^ sis was don ® t0 est imate the fluid film thickness for the 57mm 
rpflr nr O 0Per H tin * 9 n in ? x ^ 9en * The details of the analysis are contained in the 
reference and will not be repeated here. The same assumptions and reservations 

s ?f d , ir ! tho reference hold for this analysis. Figure 5.12.3 shows the fluid 
viscosity data, derived pressure coefficient, and estimated initial viscosity for 0„. 
^omnoM+ 5 * . shows ^the estimated fluid film thickness as a function of contact 
So r3t r* Uc f 01 hl . gh temperature 0„ viscosity data limited the calculation to 
, A . S was sh °'yn With N 2 , the film* thickness for 0 is too small to provide 

th^order of^r^nches ° f ^ l0aded contacts since Wc* 1 surface finishes are on 

5.12.2 Fluid Film Thickness Estimates for Material Fatigue Tester 

. Tud [J can t fil , m thickness was estimated for a material fatigue tester 

to investigate possible degradation of fatigue life due to loss of lubrication. The 
tester configuration and design characteristics are shown in Figure 5.12.5. The 
lubricant is MIL-L-7808G. The lubricant properties are also provided in Figure 
5. 12.5. The results indicate that the film thickness to roughness ratio (h/o) is not 
large enough to provide full EHD boundary lubrication. 

1 ubri rant 6 flit 1 corre r J at T° ns available in the literature for the estimation of 

minimi™ knes - s * Da ^ on M developed the following formula to calculate 

minimum film thickness in a rolling contact. 

1) 4-. 1 <w(»c') a - 6 (“? “ )0-7 

R x il'l 

Grubin developed a similar formula: 


2 ) 


■ 1-871, iV 1 ' 7 



(-rV)j 

E R x 


091 
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R 'x = 


1 


1 


jr + jr 
K 1 k 2 

2 

A = Viscosity pressure coefficient (in /lb) 
u = Ambient viscosity |1bs - sea 

\ in^ Al, + U 2 

U = Entrainment velocity = -=—* — - (in/sec) 

Wy 

W = 4 yj a = Load per unit length (lbs/in) 

Wy = Contact surface load (lbs) 
a = Semi major axis of contact ellipse (in) 

E = v Assumes properties of wheel and roller .to be equal 
E = Modulus of elasticity (lbs/in^) 
v - Poisson's ratio 


The above relationships are for isothermal conditions. H. S. Cheng (42) devel- 
oped a Grubin type inlet film thickness analysis including the thermal effects. The 
results of this work yield correction factors to the isothermal film thickness 
estimates as follows: 

3) \>in = h ' 0 s ' 0 so thermal) 


L is a reduction factor to account for the thermal effects on the lubricant and 0 

is a factor to account for side leakage in the contacts. Based on the lubricant 

properties, contact geometry, and entrainment velocity, these factors were estimated 

as follows: _ - 

0y = 0.87 


0 S = 0.9 

Dawson's film thickness formula can be adjusted to account for thermal and side 
leakage as follows: 


4) 



= 1.237 (XE* ) 6 



A similar adjustment can be made to Gurbin's formula. The minimum film thickness 
estimated from equation 4 is 3.5y inches. Gurbin's formula, including the reduction 
factors, provides a minimum film thickness of 6.2y inches. The difference is due 

primarily to the differences in exponents assigned to the load terms in the two 
equations. 


To further investigate the validity of these estimates, data from Reference (43) 
were compared with the above film thickness estimates. Figure 5.12.6 gives the 
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lubricants tested and experimental correlation factors determined from the tests. As 
shown, the type II Ester has viscosity values closest to the MIL-L-7808G lubricant. 

formula 6 fi1m thickness data fron ' refe rence 43 was correlated with the following 
h_ 


5) ^ •[<> »j l^n] 


p is a reduction factor for high stress conditions. These test data did not include 
stress conditions greater than 350ksi. Since the test fixture is stressed to 700ksi, 

woM?J U L°n cS £. u I d not be obtained directly from the test data. An optimistic value 
would be 0.5 which corresponds to a Hertz stress of 350ksi. 

Using the correlation factors for type II Ester and the MIL-L-7808G fluid prop- 
erties, the estimated film thickness from equation 5 is, 1.375^ inches. 

The surface roughness for the sample and wheels were given as: 


TEST SAMPLE 
WHEEL 


4mT0 8m INCHES 
6m INCHES 


The combined RMS surface roughness is therefore: 


CT = V^> 2 + (6u) 2 ' = lOu inch 
The results of the previous analysis are summarized below: 


CORRELATION 
Dawson 
Gurbin' s 
Eq 5 


MINIMUM FILM 
THICKNESS y INCHES 

3.5 

6.2 

1.4 


FILM THICKNESS TO 
ROUGHNESS RATIO h/a 

.35 

.62 

.14 


chnlnn 1 !! 9 t0 + .Reference (42), if h/a >4, thick EHD film exists and the lubrication 

likPlv be f s % tls ! act0 P; If Y h /f <4 > Partial EHD film exists and lubrication is most 
likely satisfactory but may deteriorate when h/a approaches unity. If h/a <1, lubriea- 

n !f r 9 1 . na1 . , and operation in this region is in general not recommended without 
nrrnr y 1 ubn catl on additives. Pitting, scoring and progressive wear are likely to 

5 * Effect of Frictio n on Rolling Element Contact and Subsurface Stresses 

affects ° f friction on the magnitude and distribution of normal and shear 
ro . 1ng element contacts has been investigated. Reference 37 contains 
linp rnntart eSS |° nS for . th ® normal and sheer stresses for rolling elements producing 
th!! ^ v C V- Lin / c ? ntact allows the stress problem to be solved in closed^ form by 
Hn tl0n t°[ i P J a r str ? ss and P 1ain strain assumptions. Although these as- 

solutions should apply f ° r P °’ nt C °" taCt ’ the trentlS 1de " t1fied in the “»f ct 
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T^ e of surface friction on the location of the maximum shear stress was 

determined by solving for shear stress and setting the derivative to zero. Figure 
b.13.1 shows the results. As shown, the location of the maximum shear stress along the 
axis in the rolling direction (x axis) is a function of the coefficient of friction. 
Also shown is the expression of the maximum surface shear (T ) as a function of 
friction. The maximum shear stress for zero friction occurs at m l x subsurface depth of 
c 1 /’ f nd stress ratio is 0.3. This point is shown on the vertical axis 
oi figure 5,13.1. As shown in the figure, the maximum shear stress moves to the 
surface for friction factors greater than 0.19. 

Shown in Figure 5.13.2 are the maximum surface principal stresses. As shown, 
these are all compressive stresses and occur in the surface as a function of the 
friction factor. The maximum tensile surface stress is shown in Figure 5.13.3, and 
occurs at (-x/a-1) in the contact surface. 


This analysis shows that the magnitude and location of the maximum contact stress- 
es can be significantly affected by the friction force in the contact surface. Since 
these stresses increase with increased friction, this effect can be especially impor- 
tant for rolling elements with marginal lubrication such as high speed bearings operat- 
ing in cryogenics. These trends may help explain the surface failures observed in 
cryogenic bearing tests as opposed to the more common subsurface failures observed in 
well lubricated bearings. 


„ This analysis has been restricted to line contact, a similar analysis for point 
contact requires evaluation of a three dimensional stress field and can best be handled 
by finite element computer codes. 

5,14 Effects of Bearing Race Curvature and Diametral C learances on Operating 
Characteristics — ~ c — a - 


** I he i" Q - X turbQ P um P turbine end bearing has been evaluated to investigate the 
effects of increased inner race curvature and increased diametral clearance on various 
operating characteristics. 

Figure 5.14.1 provides the friction heat generated per bearing pair as a function 
of radial reaction. The three curves shown represent different bearing parameters that 
were analyzed. The solid line curve represents a diametral clearance of .0036 inches, 

S:Jln a n b S nn9 Tu COnt ?f t ang - 1e of 19,95 de 9 rees • The inner and outer race curvatures 
being 0.53. The alternating long and short dashed line is for a bearing with a 
diametral clearance of .0063 inches, a contact angle of 22.89 degrees and inner and 
outer race curvatures of 0.55 and 0.53 respectively. The dashed line represents a 
bearing having a diametral clearance of .0059, a contact angle of 22.14 degrees and 
inner and outer race curvatures of 0.55 and 0.53 respectively. 

J y n i 5n was xu pe u f0rn, . ed using a shaft speed of 30 ’ 000 RPM and a friction 
coefficient of 0.20. The bearing component temperatures were assumed to be initially 

temperature of liquid oxygen. A complete thermal analysis would be 
necessary to determine the appropriate bearing component temperatures. A 1000 pound 
preload was also placed on each bearing of the pair modeled. 

Figure 5.14.2 illustrates the frictional heat generated at the inner and outer 
races per bearing pair. The frictional heat represented in Figures 5.14.1 and 5.14.2 

would nkely be more severe if the temperature distribution in the bearing were ac- 
counted for. 3 
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FIGURE 5.13.3 MAXIMUM TENSILE SURFACE STRESS vs FRICTION 
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FIGURE 5.14.1 

TOTAL FRICTIONAL HEAT GENERATED - VS - RADIAL REACTION 
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o 1000 POUND PRELOAD PER PIT 
e 30,000 RPM SHAFT SPEED 
o TEMP: -T 53* C 
® f ‘ 0.20 


7.0 







“ = 19.95' P d = .0036" 

Curvature = .53, .53 

“ * 22.89° Pd = .0063" 

Curvature - .53, .55 

u = 22.14° Pd = .0059“ 

Curvature = .53, .55 


FIGURE 5.14.2 

FRICTIONAl HEAT WiERATEO AT INNER ( OUTER RACE PER 
SEARING PAIR - VS - RAOIAl REACTION 
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The maximum stress in the contact area as a function of racha' ' rea h ctl °g,, ls 1 
c . M „ a c la o it ran he ^een that these stresses increase as the radial ioaa 

increases, as would be expected. The two cases in which the inner and ° ut ^ ra ^ 
curvature are 0.55 and 0.53 respectively appear to increase (solid line) 

coverage upon each other as the radial load increases. h 
shows a non-linear increase as a function of radial reaction. 

Fiaure 5 14 4 and 5.14.5 represent the minimum inner and outer contact angles as a 
functlJTrf 6 «d1al loading. As' shown in Figure 5.14.4, the .inner contact ang e d - 
creases nonlinearly as the radial reaction increases. Figure 5.14.5 1 1 lustrates 
? i near decrease of the outer contact angle as the radial reaction increases Also 
shown are the effects of the proposed design changes in curvature and clearances 

operating contact angles. 

The loss of internal or operating clearance as a function of radial reaction is 

radial clearance change is 64% (current clearance = 3.6 mils; modified - 5.9 mils), the 
operating clearance change is ^112%. 

aSSKT r.£ s 

stiffnesses 0 ",™ calculated to bo 1.25 x 10 6 lbs/in for the control ■ case (P. - 
nnifi'M l J153 x 10° lbs/in for the case in which the diametral clearance is .0059 and 
9 ?56 l'l0^1bs/in f7. cll.ce of .0063”. The effects of including the faring 
component temperature distribution due to thermal /rech»ni cal ' 1 "teracti on of th^bea^g 
would tend to reduce clearances and probably result in a stiffer bear _ y 
represented in Figure 5.14.7. 

5 . 15 Estimate of Cage Forces Caused by Uneven Ball Hear for Turbopump 2389R1 Turbine 
End Bearing 

An analysis was conducted to estimate the cage loads of the #3 bearing in HP0TP 
2309R1 when the webbing between the #1 and #13 PO ck . e . ts J?^ 

ranp loads were assumed to be caused by uneven ball wear. Potential inertia Torce 
caused ° by ^atin^^radial loads and ball/cage transient interactions were not con- 
sidered. The objective was to determine if quasi -steady loads 

ball wear could be severe enough to explain the cage failure. The ball wear measured 
after disassembly is shown in Figure 5.15.1. 

The larae variance in ball diameter measured (.499 in to .487 in) prevented a 
straiaht forward use of the Shaberth computer program because the equilibrium equa- 
tions 9 as coded in Shaberth, use a constant diameter ball throughout the bearing 

Therefore each ball had to be simulated separately taking into account the effect o 
thp nthpc differeit size balls in the bearing. This was done by assuming that the 
Urqest tall in ^the bearing (D = .499 in) controlled the spacing between the inner and 
iutfr racesVL ■ the bearing contains perfectly stiff races This means that the 
free contact angle based on the largest diameter ball is used (D - -499 in » ^PHA 
?3 1681 This contact angle accounts for the increased diametral clearance caused by 
wear of the largest ball 9 (twice the .001 in wear). When smaller balls are m the 
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FIGURE 5. 14.3 

MAXIMUM HERTZ STRESS - VS - RADIAL REACTION 



FIGURE S. 14.4 

MINIMUM INNER CONTACT ANGLE - VS - RADIAL REACTION 
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FIGURE 5.14.5 

MINIMUM OUIEP CONTACT ANGLE . VS RADIAL REACTION 
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FIGURE 5. 14.6 

OPERATING CLEARANCE • VS - RADIAL REACTION 
S.O, 


9 lOOO POUND PRELOAD P£R BRG 
o 30.000 RPh SKAFT SPEED 
e TO*: -1$3’C 
o f ■ 0.20 


UJ 



1.0 


* * 19.95* Pu ' .0036* 
Curvature * .53, .S3 

a « 2Z.89* P d * .0063“ 
Curvature ■ .53. .55 

-■* * 22.14* P d * .0059" 
Curvature « .53, .53 


1 1 , , 
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RADIAL REACTION PER BEARING PAIR {POUNDS) 


FIGURE 5.14.7 EFFECTS OF RACE CURVATURE A DIAMETRICAL CLEARANCES OH BEARING STIFFNESS 
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FIGURE 5.15.1 MEASURED BALL HEAR FROM FAILED BEARING 
NUMBERS IN PARENTHESES ARE AVERAGES RET1IEEN ADJACENT BALLS 


(.00575) 



TO DETERMINE THE RELATIVE RADIAL CLEARANCE FOR BALL *3 


Diametral clearance for »3 (.013" + .003" = .016'') 

- Radial Play accounted for by largest ball (2 x ,001 - ,002") 

'Relative Radial Play' associated with ball «3 



.016" 

-. 002 " 


. 014 " 


FIGURE 5.15.2 SCARING CHARACTERISTICS AS FUNCTION OF TEST BAU NUMBER 
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bearing with the largest, there is increased diametral clearance associated with each 
smaller ball. The difference between this clearance for each smaller ball and the 
clearance accounted for in the contact angle of the largest ball is considered to be 
radial^ play in the Shaberth input data. Figure 5.15.1 shows how this relative radial 
play is determined for ball #3. Figure 5.15.2 compares this clearance with the 
measured ball wear. 

Figures 5.15.3 and 5.15.4 show the data generated when each ball (simulated by its 
radial play) was run with an 850 lb radial load and an 1000 lb axial preload. Normal 
forces and orbital speeds for each azimuth location are calculated and these are taken 
to be the values for a particular worn ball at different ball positions around the 
bearing. Therefore, from these tables, the complete bearing can be constructed at 
various orientations. Marked on the tables are the values for the normal forces and 
the orbital speeds when ball #13 is directly underneath the load. The normal forces 
and orbital speeds for this orientation are plotted in Figure 5.15.2. 

Figure 5.15.5 shows the ultimate force required to break the cage at the minimum 
cross-section as a function of temperature. The tension force required to break the 
cage was taken, from this chart, to be approximately 2100 lbs. The use of this value 
assumes that, at the minimum cross-section, the cage webbing would be running at a 
relatively low temperature (-100°F) for the following reasons; 

(1) The balls contact the cage only on the pocket sides in the direction of rotation, 
not on the webbing. 

(2) The cage material is a fairly good insulator so that much of the heat generated at 
the contact would not be conducted around to the webbing, but would be carried 
away by the LGX flowing across the cage. 

The cage speed shown in Figure 5.15.2 (1904 RADS/SEC) was estimated to be that 
speed which gave the lowest net force on the cage. Since all the balls are actually 
travelling at the cage speed, the plot of the orbital speeds for each ball provides an 
indicator of which balls are pushing the cage and which are being pulled by the cage. 
The force exerted on the cage pockets is then equal to the friction force between the 
balls and the races that must be overcome if slip is to occur. This same orientation 
is shown again in Figure 5.15.6. The arrows on Figure 5.15.6 show the direction of the 
friction forces on the cage pockets. This configuration was chosen first because it 
should give the largest tension force on the pocket that failed. However, as shown, 
even if all the dragging balls were segregated from those pushing the cage, the resul- 
tant force on the cage (491 lb) would not be enough to cause the failure. 

Another orientation is shown in Figure 5.15.7, where the ball with the highest 
relative radial clearance (#10) is beneath the load. The worst cage force in this case 
is 623 lbs and this is also not enough to cause failure. 

In summa ry , the model used here does not predict a cage failure due to quasi - 
steady loads induced by uneven ball wear. It appears that there is an unstable dynamic 
interaction between the balls and the cage that is not addressed by the quasi -steady 
state analysis described herein, which produces considerably higher loads resulting in 
the cage failure. An estimate of these forces would require a detailed dynamic analy- 
sis of the cage/ball system. 
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FIGURE 5. 15.3 

SUf! OF NORMAL FORCES ON INNER AND OUTER RACE, LB 



AXIAL LOAD = 1000 LB o: = 23 168^ 

RADIAL LOAD = 850 LB d = .499" 

RESULTS FROM SHABERTH ANALYSIS OF FAILED BEARING 


FIGURE 5.15.4 

BALL ORBITAL SPEEDS, RAD / SEC 
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